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ABSTRACT 


Y 

This  report  presents  an  updated  version  of  an  earlier  simplified  mode 
conversion  program  for  VLF/LF  propagation  in  the  earth-ionosphere  waveguide. 
The  new  program  includes  the  provision  for  calculating  at  an  arbitrary  height 
within  the  guide  all  three  electric  field  components  generated  by  an  electric 
dipole  of  arbitrary  orientation  and  height  within  the  guide.  The  program  is 
designed  for  treating  air  to  air,  ground  to  air  or  air  to  ground  VLF/LF 
problems  involving  a waveguide  channel  which  is  horizontally  inhomogeneous 
along  the  direction  of  propagation. 
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I.  Introduction 

This  report  is  a continuation  of  a series  (references!,  2 and  3)  which 
describes  computer  programs  developed  for  calculating  VLF/LF  field  strengths  in 
the  earth-ionosphere  waveguide  when  allowance  must  be  made  for  horizontal 
inhomogeneity  in  the  direction  of  propagation.  Thus  these  programs  are 
particularly  relevant  to  the  problem  of  propagating  across  the  terminator  as  well 
as  to  the  problem  of  propagating  in  an  artifically  disturbed  environment.  As 
in  the  previous  works,  the  present  program  is  based  upon  a slab  model,  assumes 
waveguide  invariance  in  the  direction  normal  to  the  great  circle  path  between 
transmitter  and  receiver  and  neglects  reflections  resulting  from  inhomogeneity 
along  the  direction  of  propagation.  The  field  calculations,  principally 
through  waveguide  modal  constant  inputs,  do  allow  for  vertical  inhomogeneity 
as  well  as  anisotropy  of  the  ionosphere.  In  the  previous  works,  field  strength 
calculations  or  mode  sums  are  generated  for  the  vertical  electric  field,  Ez, 
at  the  ground  produced  by  a ground  based  vertical  dipole.  The  present  program 
differs  from  that  of  reference  3 only  to  the  extent  that  the  present  program 
can  be  used  to  calculate  all  electric  field  components  Ez,  Ex  and  E^  for  any 
receiver  height  within  the  guide  (x-z  is  the  plane  of  propagation).  Further- 
more, the  field  strength  calculations  can  be  made  for  electric  dipole  exciters 
of  arbitrary  orientation  located  at  any  height  within  the  guide.  Thus  air  to 
air,  ground  to  air  or  air  to  ground  VLF/LF  propagation  problems  involving  a 
horizontally  inhomogeneous  waveguide  channel  may  be  treated  using  the  present 
program. 

In  addition  to  familiarity  with  the  simplified  mode  conversion  model 
concepts  documented  in  references  2 and  3,  familiarity  is  assumed  with  reference 
4 which  describes  a Fortran  IV  program  for  obtaining  mode  constants  as  well  as 
excitation  factors  for  electric  dipoles  of  arbitrary  orientation  located  at 


any  height  within  the  earth  ionosphere  waveguide.  Crucial  inputs  from  the 

latter  program  are  the  ground  eigenangles,  four  independent  quantities  from 

which  a tensor  array  of  nine  excitation  factors  relating  to  end-on,  broadside 

or  vertical  dipole  excitation  of  E , E and  E may  be  determined,  and  a polari- 

z x y 

zation  factor  which  introduces  the  proper  amount  of  TE  wave  into  the  modal 
height  gains.  These  quantities  for  each  mode  and  slab  serve  as  input  for  the 
present  program. 

Principal  outputs  of  the  present  program  are  mode  conversion  coefficients 
(in  a generalized  sense)  and  mode  sum  plots  as  a function  of  distance  from 
the  transmitter  for  the  three  electric  field  components  for  four  orientations 
of  the  electric  dipole  exciter.  The  transmitter  and  receiver  must  be  within 
the  earth  curvature  dominated  portion  of  the  guide  but  otherwise  their  altitude 
is  arbitrary.  Since  the  mod  .inversion  coefficients  are  independent  of  the 
location  of  the  horizontal  inhomogeneity  relative  to  the  transmitter, 
provision  is  made  for  moving  the  inhomogeneity  in  increments  and  plotting  mode 
sums  for  the  incremented  distances  (this  option  is  useful  only  if  the  ground 
conductivity  and  geomagnetic  orientation  may  be  regarded  as  constant  over  the  path). 

In  section  II  the  mode  conversion  model  is  reviewed  and  relevant  formulas 
summarized.  A description  of  the  program  is  given  in  section  III  and  results 
presented  in  section  IV.  The  appendix  contains  a program  listing. 
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II.  Summary  of  Equations 

Inputs  to  the  present  program  for  each  slab  and  mode  are  the  ground 
eigenangles,  the  Tj  1 s defined  below  and  the  polarization  ratio  f also 
defined  below.  The  Tj's  are  readily  obtainable  from  the  waveguide  program 
documented  in  reference  4 as  are  the  eigenangles  and  polarization  ratio.  These 
quantities  are 
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S is  the  sine  of  the  eigenangle  and  F1  (en)  the  derivative  of  the  mode  equation 


evaluated  at  the  eigenangle,  ep.  The  R and  R’s  represent,  respectively,  ele- 


ments of  the  reflection  matrix  looking  into  the  ionosphere  and  towards  the 
ground  from  ground  level.  Consistent  with  the  usual  notation,  the  first  sub- 
script refers  to  the  polarization  of  the  Incident  wave  and  the  second  subscript 
refers  to  the  polarization  of  the  reflected  wave.  Equation  (5)  gives 
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the  ratio  of  to  hy  at  the  ground.  The  function  f is  also  called  F^FR  in  this 

program.  The  D-.'s  are  functions  defined  below  which  are  negated  in  the  present 
* 0 

program  by  defining 

t = D„T,  , ~t=  D22T2  , Zj  = D,xTj  (6) 

In  terms  of  the  preceding  quantities  the  excitation  tensor  elements  are 
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The  columns  relate  to  excitation  of  the  electric  field  components  Ez,  Ex  and 
E^  and  the  rows  apply  to  excitation  by  a vertical  dipole,  a horizontal  dipole 
end-on  and  a horizontal  dipole  broadside.  Recall  the  geometry  of  the  situa- 
tion is  such  that  z is  taken  positive  into  the  ionosphere,  that  positive  x is 
the  direction  of  propagation  and  that  y is  normal  to  the  plane  of  propagation. 

The  excitation  factors  must  be  supplemented  with  definitions  of  the 
height  gains.  These  along  with  the  definitions  of  the  D^'s  are 
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C = cosine  of  the  ground  eigenangle 

k = the  free  space  wavenumber 

e/ eQ  = dielectric  constant  of  the  ground 


a = the  ground  conductivity 
u = the  circular  radio  frequency 
a = the  earth's  radius 


The  functions  h-j  and  h2  are  modified  Hankel  functions  of  order  1/3  (which 
are  linearly  related  to  Airy  functions)  as  defined  by  the  Computation  Laboratory 
at  Cambridge,  Massachusetts  (reference  5)  and  the  primes  on  these  quantities 
denote  derivatives  with  respect  to  the  argument.  Equation  (20)  is  the 
modified  refractive  index  which  is  chosen  to  be  unity  at  the  ground.  The 
subscript,  o,  which  appears  on  n in  equations  (14)  and  15)  signifies  that 
Eq  (20)  is  to  be  evaluated  for  z = 0.  Similarly  the  subscript  o which  appears 
on  q in  Eq  (8)  and  Eqs.  (10)  through  (17)  signify  that  Eq.  (18)  is  to  be 
evaluated  for  z = 0.  It  should  be  pointed  out  that  f^  (apart  from  a sine  of 
the  eigenangle)  is  the  height  gain  for  the  vertical  electric  field  Ez,  f2  the 
height  gain  for  the  horizontal  electric  field  component  Ex,  and  f3  for  the 
electric  field  component  E^  which  is  normal  to  the  plane  of  propagation. 

Unlike  the  eigenangles,  the  T. 's  and  the  polarization  ratio,  the  height  gains 

J 

are  generated  in  the  present  program. 

The  final  quantities  required  for  the  mode  sum  evaluations  are  the 
generalized  mode  conversion  coefficients  which  are  also  calculated  in  the 
present  program.  Details  of  their  calculation  are  described  in  references  2 
and  3 and  we  will  remark  here  only  that  the  quantities  relate  to  the  amplitude 

X L. 

of  the  jin  mode  in  slab  p generated  by  virtue  of  a unit  amplitude  wave 
incident  in  the  transmitter  region  (slab  NTR  in  the  present  program  - see 

Fig.  1).  1 T J 


m ; m-i  i 
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Xh-1 
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XMTR 

Fig.  1 MODE  CONVERSION  MODEL 
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In  terms  of  the  excitation  factors,  height  gains  and  generalized  mode 

i,L 

conversion  coefficients  the  electric  field  components  E^p  in  the  p n slab  may 
be  written  as  follows 
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The  receiver  altitude  is  zR  and  the  transmitter  altitude  Zy.  The  final  sub- 
script on  the  x 1 s and  f's  denotesmode  indices  whereas  the  index  i takes  on 
the  values  1,  2 and  3.  Consistent  with  the  previous  definition  i = 1 Ez, 

i = 2 + E and  i = 3 + E . The  constant  Q is 

* y 

Q = 0.03248kA/F~ 

with  the  free  space  wavenumber,  k,  in  inverse  km  and  F the  frequency  in  kHz. 
The  symbol  • represents  the  Kronecker  delta.  That  is 

= 1 i = j 


The  angles  y and  <t>  determine  the  orientation  of  the  electric  dipole  source 
relative  to  the  x,  y,  z coordinate  system  as  shown  in  Fig.  2 
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Two  distinct  options  are  available  with  the  present  program.  One  is  for 
, field  calculations  (amplitude  and  phase)  as  a function  of  range  for  a fixed 

location  of  the  horizontal  inhomogeneity.  The  second  allows  for  field 
calculations  at  two  distinct  receiving  points  along  the  same  great  circle 
path  as  a function  of  position  of  the  horizontal  inhomogeneity  (this  option 
• is  useful  only  if  the  ground  conductivity  and  the  geomagnetic  parameters  are 

invariant  over  the  path).  Amplitude  is  expressed  in  dB  above  a microvolt 
per  meter  for  a one  kilowatt  radiator  and  phase  in  degrees  relative  to  free 
space.  Both  of  these  quantities  are  printed  outputs  of  the  program.  Amplitude 
plots  may  also  be  requested. 
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III.  Description  of  Program 
A.  General  Comments 


To  handle  horizontal  inhomogeneities,  the  ionosphere  is  divided  into  a 
series  of  vertical  slabs,  as  described  in  section  II.  These  slabs  are  labeled 
1,2,  ...  , M as  shown  in  Fig.  1 and  the  boundaries  between  the  slabs  have 

coordinates  X-| , y^,  ...  , X^.  For  each  slab  and  for  each  mode  the 

ground  eiqenangle,  the  T.'s  defined  by  Eqs.  (1)  through  (4)  and  the  function  f 

J 

(called  F0FR)  defined  by  Eq.  (5)  must  be  provided. 

B.  Description  of  Input 

All  input  to  the  mode  conversion  program  is  given  in  a data  deck  on  the 
standard  input  unit.  A listing  of  sample  input,  showing  the  data  deck  setup  is 
given  on  pages  17  through  ly. 

There  are  two  parts  to  the  input.  The  first  part  is  read  in  by  means  of 
a Fortran  IV  NAMELIST  input  format.  The  first  card  of  each  set  of  input  must 
contain  a blank  in  column  1 and  &DATUM  in  columns  2-7.  This  is  followed 
by  at  least  one  blank  and  then  data  items  separated  by  commas.  The  data  items 
have  the  following  forms:  (all  cards  must  have  a blank  in  column  1) 

variable  name  = constant, 

or 

array  name  = set  of  constants,  (all  separated  by  commas),  where 
successive  occurrences  of  the  same  constant  can  be  represented  by  k*  constant; 
for  example,  in  the  sample  input,  sigma  3 9*4.64  means  that  the  conductivity 
for  all  nine  slabs  is  4.64  mhos/m. 

The  following  variables  and  arrays  may  be  specified  in  the  NAMELIST 
input: 

FREQ  - frequency  in  kHz. 

SIGMA  - ground  conductivity  in  mhos  per  metre.  Note  that  a ground 

10 
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conductivity  for  each  slab  is  required  and  that  SIGMA  is  dimensioned  for  25. 

EPSL0N  - permittivity  of  ground  in  Farads  per  metre.  Note  that  a ground 
permittivity  for  each  slab  is  required  and  that  EPSL0N  is  dimensioned  for  25. 

RH0MAX  - maximum  horizontal  distance  in  km  at  which  field  strengths  are 
desi red. 

RH0MIN  - minimum  horizontal  distance  in  km  at  which  field  strengths 
are  desired. 

DELRH0  - horizontal  increment  in  km  for  which  successive  field  strengths 
are  computed. 

IDPL0T  - literal  constant  up  to  40  characters  which  is  printed  on  plots 
produced.  For  example,  in  the  NAMELIST  input  data,  place  the  card  IDPL0T  = 
'FIELD  STRENGTH  PL0T1.  If  IDPL0T  is  set  to  zero  or  if  IDPL0T  is  omitted  from 
input  no  plots  are  produced. 

NRSLAB  - Number  of  slabs  in  the  model. 

NRM0DE  - Number  of  modes  to  be  handled  in  the  program. 

NTMAX  - Number  of  times  the  transmitter  - terminator  separation  is  in- 
cremented. 

DELTAX  - Distance  in  km  by  which  transmitter  - terminator  separation  is  in 
cremented. 

XVAL  - Horizontal  position  in  km  of  boundaries  between  adjacent  slabs. 

Note  that  XVAL  can  be  negative  and  that  it  is  dimensioned  for  25. 

IFIRST  - Is  set  to  1 in  the  first  set  of  NAMELIST  input.  If  more  than 
one  set  of  input  is  used  set  IFIRST  = 0 in  the  second  set. 

LAST  - Is  set  to  1 in  the  last  set  of  NAMELIST  input.  If  the  user  has 
requested  plots  this  causes  the  end  of  file  to  be  written  on  the  plot  tape. 

IPLT0P  - Plotting  option  flag.  If  IPLT0P  = 1 two  plots  (field  amplitude 

in  dB  above  a yv/m  for  1 kw  radiated  power  versus  transmitter  - terminator 
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distance  for  two  receiver  positions)  are  produced  for  each  set  of  input.  If 
IPLT0P  = 2,  NTMAX  plots  (field  amplitude  in  dB  above  a yv/m  for  1 kw  versus 
distance  from  transmitter)  are  produced  for  each  set  of  input. 

XMIN  - minimum  value  of  X on  plot. 

XINC  - increment  of  X scale  in  km/inch. 

YMIN  - minimum  value  of  y on  plot. 

YINC  - increment  of  y scale  in  dB/inch. 

SIZEX  - size  of  X axis  in  inches. 

SIZEY  - size  of  y axis  in  inches. 

GAMMA  - dipole  orientation  angle  relative  to  z axis  (see  Fig.  2).  Note 
that  GAMMA  is  dimensioned  for  4. 

PHI  - dipole  orientation  angle  relative  to  X axis  (see  Fig.  2).  Note 
that  PHI  is  dimensioned  for  4. 

NRP  - number  of  GAMMA-PHI  pairs  up  to  4. 

ZT  - transmitter  altitude  in  km. 

ZR  - receiver  altitude  in  km. 

INTFLG  - Printing  option  flag.  INTFLG  must  be  set  to  1 if  printout  of 
height  gain  integrals  is  required.  See  reference  2 for  an  output  of  height 
gain. 

IPRNTA  - Printing  option  flag.  IPRNTA  must  be  set  for  1 if  printout  of 
generalized  mode  conversion  coefficients  is  required. 

The  end  of  the  NAMELIST  input  is  signaled  by  SEND. 

The  second  part  of  the  input  follows  the  NAMELIST  input.  Three  cards  per 
mode  are  required.  The  first  contains  the  eigenangle  at  the  ground  and  T-|  and 
T2  (see  Eqs.  1 and  2).  The  second  contains  the  eigenangles  at  the  ground 
(duplicate  input)  and  T3  and  T^  (see  Eqs.  3 and  4).  The  third  card  contains 
f (see  Eq.  5)  and  T0PHT.  The  latter  in  this  program  is  a real  variable  and 
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represents  the  height  above  ground  in  km  above  which  height  gains  are  dis- 
carded. T0PHT  is  dimensioned  for  25.  There  must  be  3*NRM0DE  cards  for  each 
slab.  The  3*NRM0DE  cards  for  slab  2 follow  those  for  slab  one  and  so  on  up 
to  slab  number  NRSLAB.  Although  ordering  of  modes  is  not  critical,  we  have 
followed  the  practice  of  ordering  them  according  to  their  real  parts.  (The 
mode  with  the  largest  real  part  is  called  mode  1). 

C.  Description  of  Output 

The  sample  output  shown  on  pages  20  through  42  begins  with  an  abbreviated 
listing  of  the  NAMELIST  input  variables.  This  is  followed  by  a printout  of 
THETA,  Tj's  and  F0FR  for  all  slabs  and  modes.  THETA  is  the  eigenangle  at  the 
ground,  the  T.'s  are  the  complex  quantities  given  by  Eqs.  (1)  through  (4)  and 

w 

F0FR  is  the  complex  quantity  f defined  by  Eq.  (5). 

The  principal  output  of  the  mode  conversion  program  begins  on  page  21 

where  the  total  mode  conversion  coefficients  defined  by  Eq.  (24)  of  reference 

2 are  printed  for  each  slab  (we  remark  that  they  are  mode  conversion  coefficients 

in  the  generalized  sense  that  the  modes  are  not  an  orthogonal  set).  The 

tabulation  represents  the  conversion  from  k to  j and  appears  as  printout 

because  IPRNTA  was  set  to  1 in  the  input.  On  the  other  hand  height  gain 

integrals  have  been  suppressed  in  the  output  because  INTFLG  was  not  set  equal 

to  1 in  the  input.  Since  IPLT0P  was  set  equal  to  2 the  sample  output  shows 

mode  sums  for  the  three  electric  field  components  E , E and  E as  a function 

z x y 

of  transmitter-receiver  distance  ranging  from  RH0MIN  to  RH0MAX  at  DELRH0  intervals. 
The  mode  sums  are  listed  in  dB/yv/m  for  a one  kilowatt  radiator  and  the  phases 
in  degrees  relative  to  free  space.  Because  NRP  = 4 in  the  input,  there  are 
four  GAMMA-PHI  pairs  (i.e.  four  antenna  orientations)  for  which  the  mode  sums 
are  computed.  Shown  on  pages  45  , 46  and  47  in  section  IV  are  reductions 
of  the  plots  generated  by  the  mode  conversion  program  for  this  case  along  with 
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WKB  results.  The  mode  conversion  results  and  WKB  results  are  shown  together 
for  program  check  purposes  as  discussed  in  section  IV. 

D.  Program  Layout 

This  subsection  describes  the  basic  features  of  the  mode  conversion 
program  listed  in  the  appendix. 

Reading  and  printing  of  input  quantities  occurs  in  MAIN  as  does  calculation 
of  constant  factors.  These  quantities  are  assigned  to  common  areas  MCINPT  or 
MCST0R.  MAIN  calls  HTINTL  and  HTGAIN  for  each  slab  MM. 

SUBROUTINE  HTINTL  (CAPI,  N0RM,  I FLG,  M,  INTFLG) 

HTINTL  calculates  the  height  integrals  defined  by  Eq.  (4)  in  reference  3. 
NORM  is  an  array  of  25  by  5 by  5 which  contains  all  combinations  of  modal 
integrals  for  the  slab  M.  Also  CAPI  is  an  array  of  25  by  5 by  5 which  con- 
tains all  combinations  of  modal  height  gain  integrals  for  the  slab  M and  the 
previous  slab  M+l.  IFLG  is  a control  flag  set  to  zero  in  MAIN  if  the  slab  M 
equals  NRSLAB.  It  is  set  to  1 if  M is  not  equal  to  NRSLAB.  INTFLG  is  a 
printing  option  flag.  It  must  be  set  to  1 if  printout  of  N0RM  and  CAPI  is 
desired.  N0RM  and  CAPI  are  assigned  to  the  common  area  CAP  in  MAIN  where  they 
are  called  TNORM  and  CAPI.  HTINTL  calls  MDHNKL 

SUBROUTINE  MDHNKL  (Z,  HI,  H2,  H1PRME,  H2PRME) 

MDHNKL  calculates  for  argument  Z two  independent  solutions  (HI  and  H2)  and 
their  derivatives  (H1PRME  and  H2PRME)  of  Stokes'  equation  by  methods  described 
in  reference  5.  MAIN  next  calls  HTGAIN. 

SUBROUTINE  HTGAIN  (Z) 

Z is  dimensioned  for  2.  Z(l)  is  set  equal  to  the  transmitter  height  ZT 
and  Z(2)  is  set  equal  to  the  receiver  height.  The  height  gain  functions  f^ , 
f£  and  f^  defined  by  Eqs.  (8),  (9)  and  (10)  respectively  are  computed  for  the 
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transmitter  and  receiver  heights.  These  are  made  available  to  MCFLD  and  MCFLD2 
through  COMMON/HTGN/. 


MAIN  next  calls  MCSTEP  for  slabs  equal  to  NTR  (transmitter  slab),  NTR-1, 

. . 1 with  allowance  for  changes  in  NTR  consistent  with  the  input  data. 

SUBROUTINE  MCSTEP  (M) 

MCSTEP  calls  for  CLINEQ  and  provides  as  its  output  the  mode  conversion 
coefficients  for  the  slabs  NTR,  NTR-1,  ...  ,1  for  all  values  of  NTR  con- 

sistent with  the  input  data.  The  mode  conversion  coefficients  defined  by 
Eq.  (24)  of  reference  2 are  printed  out  under  the  "A  = TOTAL  CONVERSION  CO- 
EFFICIENTS" label. 

SUBROUTINE  CLINEQ  (A,  B,  X,  N,  NDIM,  I FLAG,  ERR) 

CLINEQ  computes  the  solution  of  simultaneous  linear  equations  with 
complex  coefficients.  That  is  it  solves  the  matrix  equation 

A * X = B 

for  the  vector  X of  length  N,  given  the  matrix  A of  size  N by  N and  the  vector 
B of  length  N by  Crout's  L-U  decomposition  (reference  6).  The  A is  destroyed 
by  CLINEQ,  NDIM  is  an  integer  variable  which  must  be  greater  than  or  equal  to 
N.  IFLAG  is  an  integer  variable  normally  set  to  zero.  Setting  IFLAG  = 1 
bypasses  the  L-U  decomposition  of  A when  solutions  are  required  for  different 
B's.  ERR  is  a real  variable  computed  by  CLINEQ  which  indicates  the  relative 
errors  in  the  computed  solution  vector  X. 

SUBROUTINE  MCFLD 

MCFLD  called  from  MAIN  if  IPLT0P  = 1 computes  the  field  components  E^ 
defined  by  Eqs.  (22)  and  (23)  for  transmitter  height  ZT  and  receiver  height  ZR 
for  as  many  as  four  (GAMMA,  PHI)  pairs.  GAMMA  and  PHI  describe  the  orientation 
of  the  electric  dipole  source.  Calculations  are  made  for  ranges  RH0M1N  and 
RH0MAX  for  distances  between  the  transmitter  and  the  start  of  the  horizontal 

15 


inhomogeneity  ranging  between  XVAL  ( NRSLAB- 1 ) and  NTMAX*DELTAX  +XVAL(NRSLAB-1 ) 
at  intervals  of  DELTAX  using  Eqs.  (22)  and  (23).  Field  amplitude  outputs  are 
in  dB  above  a pv  per  metre  for  1 kw  radiated  power  and  phase  angles  are  in 
degrees  relative  to  free  space  phase. 

SUBROUTINE  MCFLD2 

MCFLD2  called  from  WIN  if  IPLT0P  = 2 computes  the  field  components  E^* 
defined  by  Eqs.  (22)  and  (23)  for  transmitter  height  ZT  and  receiver  height  ZR 
for  as  many  as  four  (GAMMA,  PHI)  pairs.  Calculations  are  made  for  transmitter- 
receiver  distances  ranging  from  RH0MIN  to  RH0MAX  at  DELRH0  intervals  using 
Eqs.  (22)  and  (23)  for  a fixed  horizontal  inhomogeneity.  Field  amplitude 
outputs  are  dB  above  a pv  per  metre  for  1 kw  radiated  power  and  phase  angles 
are  in  degrees  relative  to  free  space  phase. 

SUBROUTINE  MCPLTS 

MCPLTS  generates  six  plots  (three  field  component  amplitudes  in  dB  above  a 
pv  per  metre  for  1 kw  radiated  power  versus  distance  between  transmitter  and 
start  of  the  horizontal  inhomogeneity  for  two  receiver  ranges).  As  many  as 
four  (GAMMA,  PHI)  pairs  are  possible  so  that  each  plot  can  contain  as  many  as 
four  curves. 

SUBROUTINE  MCPLT2 

MCPLT2  generates  three  plots  (three  field  component  amplitudes  in  dB  above 
a pv  per  metre  for  1 kw  radiated  power  versus  transmitter  receiver  distance 
for  a single  location  of  the  horizontal  inhomogeneity).  As  many  as  four  (GAMMA, 
PHI)  pairs  are  possible  so  that  each  plot  can  contain  as  many  as  four  curves. 

SUBROUTINE  MAGANG  (ARG,  MAG,  ANGLE) 

MAGANG  converts  complex  number  ARG  to  polar  form  with  ANGLE  in  degrees. 
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SAMPLE  INPUT 


GDATUM 

IDPLOT=*NP« IME  = 70  TO  NPR1ME  = 71 
FREQ=21.794, 

RHOMIN=25. ,RHOMAX=50JO. ,0ELRH0=25. ,0ElTAX=0.,NTMAX=l, 

NRSLAB=9»  NRMClOE  = 5 > 

SIGMA  = 9*4. 69,  EPSLUN  = 9*. 71 72014E-09, 

XVAL=1025., 837.5,  71 2 . 5 , 58  7. 5, 4o2.  5 , 33  7.  a,  2 12  . 5 , 2b  . , 0.  , 

GAMMA  * 0.  ,90.  ,90.  , 45  • , PH  1 = 0 • , 0. , 90. , 45 . , NR P=4, 

ZT=5.,2R*  10., 

XMIN=C.tXl NC=500. «4'MIN  = -80.,VINC  = 20.,SIZEX  = 10.*SIZEV  = 8. , 

I PRNT  A= 1 , 1PLT0P=2, IF IRST*1 , 

LAST=l, 

GEND 

1 89.43093  -3.45987  4 . 329734300-03-2 . 11 976o36D-02- 1 . 41204 52  ID-1 4-2 . 46d472660-l4 

2 89.43090  -3.45987  2.175780870-08  l . 890 l a 1640-08  1.310387010-01  7. 290033340-01 

89.43090  -3.45987  4.612300-02  1.77108  -o. 547 1 1 8400-07  l . 16015 319D-0o  71.00 

1 87.19352  -1.77741  4. 8 25 1 5 7 780-05- l. 8420 2 51 80-05-6. 04653 5 120-1 1 -2 . 40 7 738400- l l 

2 87.19352  -1.77741  -4 . 263 745280-08- 5. 2 1 2643 140-08  l. 36581 734D-01  7.285013610-01 

87.19352  -1.77741  1.09558D-04  2.77337  -♦. I 1288203D-04  -1. 2373221 10-33  71.00 

1 82.32644  -1.03413  1.511146650-03-3.752447000-02  2 . 8222667 70-14-3. 3 7 2784000-1 3 

2 82.32644  -1.03410  7.596839730-08  1.068389410-07  1.547947470-01  7.2 J226265D-01 

82.32644  -1.03410  7.83611D-02  1.60618  -2.761173930-06  2.135698020-06  71.00 

1 79.02547  -1.34370  l . 1 786685 10- 04- 1 . 03 38 1 56 7D-05- 3. 26606 199D- 1 0- 9. 72 733909D-1 l 

2 79.02547  -1.34370  - 1 . 2 51 1 401 50-07- 1 . 9672 57380-0 7 1 . 7462 0 1 82D-01  7. 2 1086723D-0 1 

79.02547  -1.34370  2.424100-04  3.04333  -9. 05025258D-04  - l . 74997603D-03  71.00 

1 75.47092  -1.35774  l . 5 827 16243-03- 3 . 38o069990-02  2 . 468936600-1 3-3. 541 508280-12 

2 75.47092  -1.35774  2. 225 1 36400- 07  3. 343343850-07  2. 05880 77 ID-0 1 7.180996130-01 

75.47092  -1.35774  6.751710-02  1.60522  -9.536959350-06  7 . 0 1 72 l 882D-06  ,71.00 

1 89.42754  -3.44563  4. 3 423833 1 3-03-2 . 13 1 7 1 1 800- 02- 1 . 3322 746 10-1 4-2. 4038 5 1 230-1 4 

2 89.42754  -3.44563  2.  1 19046920-08  1.8663394  70-08  1.330  342780-31  7. 3 805 7 1 94D-0 1 

89.42754-3.44563  ■♦.63765D-02  1.77055  -e. 4o0 777220-0 7 1.125667470-06  70.88 

1 87.16147  -1.76119  4. 695 543950-03- l . 7325 91240-05 -6. 1 03696770- l 1-2. 423655940-1 1 

2 87.1614/  -1.  761  19  -4.  16237 489D  — Od—  5.  1451 5 a64D- 08  1.3  852  7 7180-01  7.3  766  07840-01 

87.16147  -1.76119  1.061610-04  2.78506  -4. 24343247D-04  -l . 2 323451 23-0 3 70.88 

1 82.29943  -1.03634  1 . 50 1040260-03-3 . 7592 3 7960-02  3. 1 18308740- 1 4-3. 3005 1333D- l 3 

2 82.29943  -1.03634  7. 433049o9D-0d  1.054259470-07  1.551452550-01  7. 33335067D-0 1 

82.29943  -1.03634  7.851200-02  1.60581  -2 . 721 160080-06  2.085929d80-06  70.88 

1 78.99378  -1.34503  1 . 1 491 903 5D-04- 9. 2202 4502 0-06-3. 28887 1380-1 0-9. 769304600- 1 1 

2 78.99378  -1.34533  - l . 2 22 1994 10-0 7- 1 . 9406608 ID-0 7 1.  773261630-31  7. 30 777892D-01 

78.99378  -1.34503  2.361150-04  3.05240  -9. 22 1 101 71 0-04  -1.762719940-33  70.88 

1 75.43374  -1.3o379  l . 603 lb53 3D- 03- 3 . 3920 35650-02  2 . 8667094 ID-l 3-3 . 465999720- 12 

2 75.43374  -1.3o379  2. 169576250-07  3.299438840-07  2.076356480-01  7.263593410-01 

75.43374  -1.36379  6.7ol53U-02  1.60566  -9. 403 71 1090-06  6.840538870-06  70.88 

1 89.42413  -3.43135  4.  3 55064960-03- 2 . 143658070- 02 - 1 . 2565007 50 -14-2. 3 39916500-14 

2 89.42413  -3. •♦3135  2 . 0638392  50-08  1 . 842 1 8t>4 70- 08  1 . 3 51 60 0070 -01  7.470059  570-01 

89.42413  -3.43135  4.663000-02  1.77003  -O.374594360-07  1 . 0922 7 lo 70-06  70.75 

1 87.12937  -1.7*541  4. 5c7999340 -05 - l . 6293 9 3290-05- 6.  161 136910- 1 1-2. 44 37433 70-1 1 

2 87.12937  -1.74541  -4 . 063 5bl990-08- 5. 0785 78300-08  1.406111950-01  7.467098770-01 

87.12937  -1.74541  1.028610-04  2. 79597  -4. 3 74240600-04  -1.267804030-03  70.75 

1 82.27203  -1.03933  1 . 5 13964290-03- 3 . 76  37 84720-02  3. 388 1 1 91 50- 1 4- 3. 23028 5920- 1 3 

2 82.27200  -1.03930  7. 259834420-06  1 . 04 l 3 5o60D-0 7 l . 5 712895 70-0 1 7. 42 7 766200-01 

82.27200  -1.03930  7.859790-02  1.60608  -2. 664848640-06  2.03686222D-06  70.75 

1 78.96090  -1.34548  l . 120422670-04- 7. 6936 8o490- 06- 3. 3 l J285440- 1 0- 9. 847235680- 1 1 

2 78.96396  - 1 . 34548  - 1 . 1 930652 50-0 7- l . 91 a6 58280-0 7 1.787470910-01  7. 4033 1 444D-0 1 

78.96096  -1.34548  2.299560-04  3.06387  -9. 429222540-04  -1.775404610-03  70.75 

1 75.39646  -1.36988  1 . 62360 7820-03- 3 . 3980 30050-02  3 . 2 3724 8290-1 3- 3. 39 1 881 29D- 12 

2 75.39648  -1.36988  2.  1 1550838D-0  7 3 . 2595 0655D-07  2 . 095 72 7700-0 1 7.384866790-01 

75.39648  -1.36988  6.771390-02  1.60608  -9.273700980-06  6.668796000-06  70.75 

1 89.42067  -3.41703  4. 36780 l 590-03-2 . 1555 97040-02- l. 18448522D-1 4-2. 276774930- 14 

2 89.42067  -3.41703  2. 0 1010536D-08  1.818570540-08  1.374169710-01  7.558491620-01 
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89.4206/  -3.41703 

1 87.0972b  -1.73005 

2 87.09728  -1.73005 
87.09728  -1. 73005 

1 82.24464  -1.04249 

2 82.24464  -1.04249 
82.24464  -1.04249 

1 7d. 92935  -1.34675 

2 78.92935  -1.34675 
78.92935  -1.34675 

1 75.35909  -1.37587 

2 75.35909  -1.37587 
75.  35909  - 1 . 375b  7 

1 89.41716  -3.40266 

2 69.41716  -3.40266 
89.41716  -3.40266 

1 87.06508  -1.71514 

2 87.06508  -1.71514 
87.06508  -1.71514 

1 82.21716  -1.04550 

2 82.21716  -1.04550 
82.21716  -1.04550 

1 78.89694  -1.34785 

2 78.89694  -1.34785 
78.89694  -1.34785 

1 75.32087  -1.38223 

2 75.32087  -1.38223 
75.32087  -1.38223 

1 89.41454  -3.39001 

2 89.41454  -3.39001 
89.41454  -3.39001 

1 87.03294  -1.70078 

2 87.03294  -1.70078 
87.03294  -1.70078 

1 82.18966  -1.04853 

2 82.18966  - 1.04853 
82.18966  -1.04853 

1 78.86419  -1.34868 

2 78 . d64 19  -1.34868 
78.86419  -1.34868 

1 75.28492  -1.38792 

2 75.28492  -1.38792 
75.28492  -1.38792 

1 89.41093  -3.37555 

2 89.41093  -3.37555 
89.41093  -3.37555 

1 87.00067  -1.68670 

2 87.00067  -1.68670 
87.00067  -1.68670 

1 82.16221  -1.05178 

2 82.16221  -1.05178 
82.16221  -1.05178 

1 78.83179  -1.34989 

2 78.83179  -1.34939 
78.83179  -1.34989 

1 75.24741  -1.39414 

2 75.24741  -1.39414 
75.24741  -1.39414 

1 89.40604  -3.35923 

2 89.40604  -3.35923 
89.40604  -3.35923 


4.688340-02  1.76951  -6.288805040-07  1.059932750-06  70.63 

4.441 751000-03-1.5291 70 71 0-05 -6. 2 14900250 *1 l -2 .462 706 130- l 1 
-3.966465770-08-5.016873020-08  1 . 43 1 122280-0 l 7. 548 702270-01 
9.96260D-05  2.80697  -4. 307255180-04  -l . 284653220-03  70.63 

1.527754700-03-3.768758320-02  3.646065450-14-3.159154480-13 
7.0b7286970-0b  l . 0292 88490-07  1 . 596239940-01  7.512977580-01 
7.869270-02  1.60636  -2. 6 5052 02 1 0-06  l . 9d 798 140D-06  70.63 

1.091984290-04-6.538902380-06-3. 336740410-10-9.890797890-11 
-1. 165503390-07-1.890994420-07  1.817086530-31  7.498799/6D-01 
2.239450-04  3.07259  -9.601867100-04  -1 . 78920 169D-03  70.63 

1.643696850-03-3.403994490-02  3. 53821 78 70-1 3- 3. 31 4067460-12 
2.063571590-07  3.216  7 90430-07  2.  1 1491 9050-0 1 7.483378540-01 
6.781170-02  1.60650  -9.136016300-06  6.503359530-06  70.63 

4.3  80692620-03-2.  16 75 3 J03O-02 - 1 . 1 160 1 1410-1 4-2 . 2 1 459 7 130-1 4 
1.957808150-08  1 . 7952 32 lbD-08  1.398035580-01  7. 64584657D-01 
4.713680-02  1.76901  -6.203482710-07  1.028617130-06  70.50 

4.31868966D-05-1.43578363D-05-6.277947520-11-2.481362760-11 
-3. 872261380-08-4.952447840-08  1.454045030-01  7.636470610-01 
9.651200-05  2.81756  -4. 64087034D-04  - 1. 301 03 72 ID-03  70.50 

l. 5 40 61 88 7D -03-3. 77331583 0-02  3.881863180-14-3.090937710-13 
6.922065980-08  1.016849470-07  1.618444330-01  7.604610030-01 
7.877860-02  l.60o62  -2. 615582 500-06  1 .94127072D-06  70.50 

1.064363780-04-5.315815950-06-3. 360 740090-1 0-9. 952781380- l l 
-l. 138402000-07-1. 86703742D-07  l . 8392 1 0800-01  7. 59 182693D-0 1 
2.181140-04  3.08246  -9.795093260-04  -1.80305 60 9D-03  70.50 

1.666511590-03-3.410825100-02  3. 83897 7980 -1 3- 3. 24 5080470-12 
2.0 1 61 4036D-07  3. 1 8 13 3 7180- 07  2. 1 32297380-01  7.571759030-01 
6. 79263D-02  1.60698-9.016843580-06  6.351563150-06  70.50 

4.3  8 7896040-03-2. 1 765 91 060-02- 1 . 0468 746 70- 1 4-2. 14  7844920-14 
1.902140900-08  1 . 7729 991 8D-08  1.438232790-01  7. 71 326704D-0 1 
4.732780-02  1.76852  -6.134690450-07  9.975804700-07  70.38 

4.201166820-05-1.347657930-05-6.338101980-11-2.500689950-11 
-3.781787040-08-4.894907260-08  1.482274170-01  7.715752830-01 
9.355630-05  2.82811  -4.773041770-04  - 1 . 3 1823 6260-03  70.38 
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TOTAL 

CONVERSION 

COEFFIC IENTS 

SLAB  NUMBER  = 

J 

* 

1 

K 

X 

l 

A= 

0. 3202  5D  00 

0.668610  00 

J 

= 

1 

K 

= 

2 

A = 

-0.120220-01 

-0.505130-01 

J 

= 

l 

K 

■ 

3 

A = 

0.19543D-02 

-0. 570360-03 

J 

s 

1 

K 

= 

4 

A = 

-0.621450-03 

0.202400-02 

J 

* 

1 

K 

= 

5 

A = 

0 . 45805  D-03 

0.  18273D-03 

J 

= 

2 

K 

X 

I 

A = 

-0. 102610-03 

0.  116010-03 

J 

= 

2 

K 

X 

2 

A = 

-0.341100  00 

0.415040  00 

J 

* 

2 

K 

3 

3 

A = 

0. 156380-04 

-0.258600-04 

J 

= 

2 

K 

X 

4 

A = 

0.81234D-03 

-0.410840-03 

i j 

X 

2 

K 

X 

5 

A = 

0.619370-05 

0.  167650-05 

j 

= 

3 

K 

X 

I 

A = 

- 0 . 39335D-02 

0.197570-02 

j 

= 

3 

K 

X 

2 

A = 

0.588830-02 

-0.783210-03 

j 

X 

3 

K. 

= 

3 

A= 

0.310320  00 

-0.790640-01 

j 

= 

3 

K 

* 

4 

A = 

-0.249970-02 

-0.856210-02 

j 

3 

K 

= 

5 

A» 

-0.628460-03 

-0. 114470-03 

j 

* 

4 

K 

* 

I 

A = 

0.243760-05 

-0. 165480-05 

j 

= 

4 

K. 

= 

2 

A = 

-0.250790-02 

-0. 558980-03 

- 

j 

= 

4 

K 

S 

3 

A = 

-0.919300-05 

0.  135260-04 

j 

= 

4 

K 

X 

4 

A* 

-0 . 87304D-0 1 

-0.100690  00 

j 

* 

4 

K 

X 

5 

A* 

0.610990-05 

0.902520-06 

j 

* 

5 

K 

X 

1 

A = 

- 0. 8 1859 D- 03 

—0. 8j6060-03 

j 

= 

5 

K 

X 

2 

A = 

0.209840-02 

0. 1409o0-02 

; 

j 

* 

5 

K 

X 

3 

A* 

0.66462D-03 

0.373030-03 

j 

x 

5 

K 

X 

4 

A* 

0.121520-02 

0.460090-03 

i 

j 

X 

5 

K 

X 

5 

A = 

-0 • 34434D-0 l 

-0.496220-01 

ii 

X 

TOTAL 

CONVERSION 

COEFFIC IENTS 

SLAB  NUMBER  = 

I 

j 

X 

1 

K 

X 

I 

A* 

0.392100  00 

-0.677470  00 

j 

X 

1 

K 

X 

2 

A* 

-0.31125D-01 

0.388720-01 

j 

X 

1 

K 

X 

3 

A* 

-0.214330-02 

-0.  176770-02 
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J 3 
J = 
J 3 
J 3 
J 3 
J = 
J = 
J 3 
J = 
J = 
J 3 
J = 
J 3 
J 3 
J * 
J = 
J = 
J 3 
J 3 
J 3 
J = 
J 3 
A = 
J = 
J 3 
J 3 
J 3 
J * 
J * 
J = 
J * 


1 

K 

3 

4 

A3 

O. 221 1 9D-02 

-0. 116350-02 

1 

K 

= 

b 

A = 

-0.107620-03 

-0.548390-03 

2 

K 

* 

i 

A = 

0.14835D-03 

0.359350-05 

2 

K 

S 

2 

A = 

3.601930  00 

-0. 731100-01 

2 

K 

* 

3 

A = 

-0.3403 10-04 

0.  lb6040-04 

2 

K 

3 

4 

A = 

-3.847860-03 

0. 250253-04 

2 

K 

= 

5 

A = 

- 3. 51 536D- 05 

-0.561920-05 

3 

K 

3 

1 

A = 

0.690950-02 

0.  1565*0-02 

3 

K 

= 

2 

A3 

-3.788310-02 

-0.263923-02 

3 

K 

* 

3 

A = 

-0.387600  00 

0.  731883-01 

3 

K 

= 

4 

A = 

0. 354520-02 

0.  1*1063-0 1 

3 

K 

= 

5 

A = 

0.925730-03 

3.406570-03 

4 

K 

= 

i 

A = 

-0.502740-05 

-0.293003-05 

4 

K 

= 

2 

A3 

0.253930-02 

0.3*5053-0 2 

4 

K 

3 

3 

A3 

0.101820-04 

-0. 2*52b3-04 

4 

K 

3 

4 

A 3 

0.192700  00 

0.254220-01 

4 

K 

= 

5 

A3 

-0.773620-05 

-0.39U13J-05 

5 

K 

= 

1 

A = 

- 0 . 14762D-02 

0.  864290-03 

5 

K 

3 

2 

A 3 

0.155950-02 

-0.377143-02 

5 

K 

= 

3 

A3 

-0.151210-02 

-0.  U4300-02 

b 

K 

3 

4 

A3 

-0.293660-02 

3.  100330-02 

b 

K 

3 

b 

A = 

0.592960-01 

-0. d l 7 100-0 1 

TOTAL 

CONVERSION 

COEFUC  I ENTS  SLAB  NUMBER  3 

1 

K 

3 

1 

A* 

3.74853D  00 

-0.316010  00 

1 

K 

3 

2 

A = 

-0.434430-01 

0.  1381*0-01 

X 

K 

r 

3 

A* 

-0.13717  D-02 

-0.311170-02 

1 

K 

S 

4 

A* 

0.185030-02 

-0. 133210-03 

1 

K 

3 

5 

A* 

0. 17119 D- 03 

-0.443360-03 

2 

K 

3 

1 

A* 

0.107740-03 

0.777620-04 

2 

K 

3 

2 

A * 

3. 6 06 94 D 00 

0.247230  00 

2 

K 

3 

3 

A* 

-0.40345D-04 

0.674430-05 
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J 

X 

2 

K = 

4 

A* 

-0.316410-03 

-0.330030-03 

J 

X 

2 

K * 

5 

A* 

-0.162890-06 

-0.612650-05 

! j 

X 

3 

K *= 

1 

A3 

0.59591D-02 

0.578120-02 

j 

X 

3 

K * 

2 

A = 

-0. 72836D-02 

-0.564410-02 

j 

X 

3 

K = 

3 

A* 

-0.44692D  00 

0.669860-01 

j 

X 

3 

K * 

4 

A- 

0.634190-02 

0.141250-01 

j 

s 

3 

K * 

5 

A = 

0.12215D-02 

0.  117200-03 

J 

X 

4 

K * 

1 

A* 

-0.301260-05 

-0.561400-05 

X 

4 

K «= 

2 

A* 

0.284200-04 

0.415570-02 

j 

X 

4 

K * 

3 

A* 

0.104710-04 

-0.271990-04 

j 

X 

4 

K = 

4 

A* 

0.233990  00 

-0.867150-01 

j 

X 

4 

K x 

5 

A* 

-0.234640-04 

0. 183640-05 

j 

X 

5 

K » 

i 

A* 

-0.176610-02 

0.830280-04 

* j 

X 

5 

K x 

2 

A* 

0.300370-02 

-0.279150-02 

; j 

X 

5 

K x 

3 

A* 

-0. 16473D-02 

-0.156700-02 

j 

X 

5 

K x 

4 

A* 

-0.375260-02 

0.324860-02 

j 

X 

5 

K x 

5 

A* 

-0.86601D-01 

-0.113360  00 

A 

X 

TOTAL 

CONVE 

RS  ION 

COEFFICIENTS  SLAB  NUMBER  * 4 

J 

X 

1 

K x 

1 

A* 

0.817960  00 

0.205590  00 

i J 

X 

i 

K = 

2 

A* 

-0.390130-01 

-0. 108330-01 

j 

X 

1 

K x 

3 

A* 

0 . 60306  0-06 

-0.365680-02 

J 

X 

1 

K * 

4 

A = 

0.630000-03 

0.  739180-03 

j 

X 

1 

K = 

5 

A* 

0.44747D-03 

-0.  U018D-03 

j 

X 

2 

K = 

i 

A* 

0.39390D-04 

0.108450-03 

j 

X 

2 

K = 

2 

A* 

0.444420  00 

0.550940  00 

j 

X 

2 

K * 

3 

A* 

-0.413720-04 

-0.473950-05 

i j 

X 

2 

K « 

4 

A* 

0.445250-03 

-0.626480-03 

j 

X 

2 

K = 

5 

A* 

0.57154D-05 

-0.730860-05 

j 

X 

3 

K = 

1 

A* 

0.287370-02 

0.879550-02 

fi  j 

X 

3 

K x 

2 

A* 

-0.477640-02 

-0.833040-02 
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J 

S 

3 

k 

* 

3 

A- 

-J. 519400  00 

0.591320-01 

J 

S 

3 

k 

X 

4 

A = 

0.987 14D-02 

0.  153110-01 

J 

= 

3 

k 

X 

5 

A* 

0.93620D-03 

-0.283090-03 

J 

= 

4 

k 

= 

1 

A = 

•0.441 16(3-08 

-0. 576080-05 

J 

= 

4 

k 

= 

2 

A* 

-0.294730-02 

0.427090-02 

J 

* 

4 

k 

x 

3 

A = 

0.90101 D-05 

-0.317060-04 

J 

= 

4 

k 

* 

4 

A = 

0.211950  00 

-0.240780  00 

J 

= 

4 

k 

* 

5 

A = 

-0.281930-04 

0.332830-04 

J 

= 

5 

k 

* 

1 

A = 

-0. 122580-02 

-0.723040-03 

J 

* 

5 

k 

s 

2 

A * 

0.408570-02 

-0.268410-03 

J 

= 

5 

k 

= 

3 

A * 

-0.187650-02 

-0. 180590-02 

J 

= 

5 

k 

= 

4 

A- 

-0.1 7384D-02 

0 . 7^6 160-02 

J 

= 

5 

k 

s 

5 

A * 

-0.188290  00 

0. 732630-01 

A 

= 

TOTAL 

CJNVERS I ON 

COEFFICIENTS  ! 

SLAB  NUMBER  = 

J 

= 

1 

k 

x 

1 

A* 

0.550330  00 

0.661000  00 

J 

« 

1 

k 

= 

2 

A = 

- J. 227360-01 

-3.259730-01 

J 

l 

k 

= 

3 

A * 

0. 147920-02 

-0  * 368 l 50-02 

J 

= 

1 

k 

= 

4 

A* 

-0 . 12503  D-Q  2 

0. 13454 J-U2 

J 

= 

1 

k 

= 

5 

A = 

J. 625560-03 

J.  150880-03 

J 

= 

c 

k 

* 

1 

A = 

•0.234090-04 

0.931160-04 

J 

= 

2 

k 

x 

2 

A = 

0.133630  00 

0 . 7 33  3o0  00 

J 

= 

2 

k 

* 

3 

A = 

-0.371450-04 

- J. 148780-04 

J 

= 

2 

k 

* 

4 

A = 

0.127050-02 

-0.918450-03 

J 

= 

2 

k 

= 

5 

A- 

0.  /81  750-05 

-0.807340-05 

J 

= 

3 

k 

= 

1 

A = 

-3.133100-02 

0. 934900-02 

J 

* 

3 

K 

= 

2 

A- 

-0.  ld786l)-02 

-0.940590-08 

J 

= 

3 

k 

X 

3 

A = 

-0.600690  00 

0.494270-0 1 

J 

= 

3 

k 

= 

4 

A = 

J.  136830-01 

0.  148740-01 

J 

= 

3 

k 

X 

5 

A = 

0.292200-03 

0. 148670-03 

J 

x 

4 

k 

= 

1 

A* 

0.254360-05 

-0.357080-03 

J 

3 

4 

k 

x 

2 

A- 

-0. 574750-08 

0.  311740-02 
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J = 4 K = 3 A*  0.53667D-05  -0.345500-04 
J = 4 K * 4 A*  0.910690-01  -0.402510  00 
J = 4 K » 5 A=  0.976610-05  0.643780-04 
J x 5 K = 1 A*  0.154830-03  -0.170990-02 
J x K = 2 A = 0.512660-02  0.202250-02 
J = 5 K * 3 A=  -0.191220-02  -0.218010-02 
J = 6 K = 4 A*  0.451 4 1 D-02  0.914880-02 
J = 5 K x 5 A = 0.283670-01  0.265570  00 

A x total  conversion  coefficients  slab  NUMBER  * 6 

J » i K = 1 Ax  0.323460-01  0.906070  00 
J = 1 K = 2 A*  -0.411230-02  -0.273840-01 
J x i K = 3 Ax  0. 2461 5D-02  -0.278110-02 
J = 1 K = 4 A*  -0.314200-02  0.151150-02 
J x i K = 5 A»  0.458340-03  0.238190-03 
J = 2 K x i A=  -0.572620-04  0.467450-04 
J x 2 K = 2 A=  -0.264140  00  0.763630  00 
J = 2 K = 3 A=  -0.285660-04  -0.210470-04 
J = 2 K = 4 A=  3.19570D-02  -0.118360-02 
J = 2 K = 5 A=  0.490450-05  -0.344850-05 
J = 3 K = 1 A*  -0.454380-02  0.7b9830-0 2 
J = 3 K = 2 A=  0.116890-02  -0.841630-02 
J x 3 K = 3 A=  -0.694490  00  0.379620-01 
J = 3 K = 4 A=  0.164460-01  0.106600-01 
J x 3 K = 5 A = 0.610550-03  0.131660-02 
J = 4 K = 1 A=  0.493000-05  0.667090-06 
J = 4 K * 2 A*  - 0. 763 19 D-02  0.136970-02 
J x 4 K = 3 A=  0. 52563D-06  -0. 335050-04 
J x 4 K = 4 A=  -0.1 482 7 D 00  -0.51006J  00 
J = 4 K * 5 A*  0.7184&0-04  J.39197D-04 
J = 5 K = 1 Ax  0.587850-03  -0.226680-02 


J 

* 

5 

K 

* 

2 

A = 

0.436340-02 

0. 309870-0 2 

J 

■ 

5 

K 

= 

3 

A = 

-0.229080-02 

-0.273560-02 

J 

3 

5 

K 

* 

4 

A = 

0 • 112  51 D-  01 

0.488890-02 

J 

= 

5 

K 

* 

b 

A= 

0.40306D  00 

0.667000-01 

A 

* 

tOTAL 

CONVERSION 

COEFFICIENTS  ! 

SLAB  NUMBER  = 

J 

= 

I 

K 

* 

1 

A = 

-0.538580  00 

0.77454D  00 

J 

= 

I 

K 

S 

2 

A = 

0.80052D-02 

-0. 179770-01 

J 

= 

1 

K 

* 

3 

A = 

0 .2  5998  D- 02 

-0. 1526*0-02 

J 

= 

1 

K 

= 

4 

A = 

-0.41 63 3D- 02 

0. 14959 J-02 

J 

Z 

l 

K 

= 

5 

A = 

-0.121690-03 

0.328310-03 

J 

= 

2 

K 

s 

1 

A = 

-0 . 5193  3 D-04 

0.505900-05 

J 

= 

2 

K 

* 

2 

A = 

- 0.6 5604 D 00 

0.607390  00 

J 

r 

2 

K 

= 

3 

A= 

-0.177610-04 

-0.209880-04 

J 

= 

2 

K 

= 

4 

A = 

0.211340-02 

-0. 134540-02 

J 

= 

2 

K 

= 

5 

A = 

0.389280-05 

0.  744470-05 

J 

s 

3 

K. 

= 

l 

A= 

-0. 58830D-02 

0.432300-02 

J 

= 

3 

K 

= 

2 

A = 

0.30027  D-02 

-0.611 110-02 

J 

* 

3 

K 

* 

3 

A= 

-0.803170  00 

0.246010-01 

J 

= 

3 

K 

= 

4 

A = 

0.162540-01 

0. 54558D-02 

J 

= 

3 

K 

= 

b 

A = 

0.237 19D-02 

0.  148430-02 

J 

4 

K 

s 

1 

A = 

0.443810-05 

0.  420020-05 

J 

= 

4 

K 

r 

2 

A = 

-0.747830-02 

-0.850940-03 

J 

= 

4 

K 

z 

3 

A = 

-0.32726D-05 

-0.286770-04 

J 

* 

4 

K 

* 

4 

A = 

-0.488750  00 

-0.479460  00 

J 

= 

4 

K 

b 

A = 

0.78476D-04 

-0. 392660-04 

J 

* 

5 

K 

= 

I 

A = 

0.739950-03 

-0.962090-03 

J 

* 

5 

K 

z 

2 

A = 

0.103150-02 

0. 363160-02 

J 

= 

5 

K 

= 

3 

A= 

-0.263180-02 

-0.320570-02 

J 

* 

5 

K 

- 

4 

A = 

0. 122060-01 

-0.408860-02 

J 

s 

5 

K 

s 

5 

A> 

0.23669D  00 

-0.533080  00 

A = TOTAL  CONVERSION  COEFFICIENTS  SLAB  NUMBER  = 8 


J 

X 

1 

K 

X 

1 

A* 

0.99/ OB D 00 

0. 20b  320-0  3 

J 

X 

1 

K 

X 

2 

A = 

-0.1OA27D-OI 

0. 31A820-02 

J 

s 

1 

K 

* 

3 

A* 

-0. 19052D-02 

0. 199060-0A 

J 

* 

l 

K 

X 

4 

A = 

0.A27380-02 

-0. 1A5220-03 

J 

X 

1 

K 

* 

5 

A = 

0.83337D-03 

-0.363280-03 

J 

X 

2 

K 

X 

1 

A = 

0 .209 1 30-0 A 

0.  1885AJ-0A 

J 

X 

2 

K 

= 

2 

A = 

3. 101 02 D 01 

-0.A682A0-02 

J 

X 

2 

K 

= 

3 

A = 

0.A9692D-05 

0. 162880-0A 

J 

X 

2 

K 

= 

A 

A = 

-0. 198910-02 

0. 998070-03 

J 

X 

2 

K 

X 

5 

A = 

-0.639060-05 

-0. 17o560-0A 

J 

* 

3 

K 

X 

I 

A = 

0.AA515D-02 

-0.265860-03 

J 

X 

3 

K 

X 

2 

A* 

-0.321030-02 

0.256260-02 

J 

X 

3 

K 

= 

3 

A* 

0.998AAD  00 

-0.556770-03 

J 

X 

3 

K 

X 

A 

As 

-0.130080-01 

0.AA7650-03 

J 

X 

J 

K 

X 

5 

A* 

-0.361720-02 

0.733230-03 

J 

* 

A 

K 

X 

I 

A* 

-0.360670-05 

-0. A163AO-05 

J 

= 

A 

K 

X. 

2 

A* 

0.666980-02 

0.  160700-02 

J 

* 

A 

K 

X 

3 

A* 

0.893A60-05 

0. lb A980-0 A 

J 

X 

A 

K 

X 

A 

As 

0.100680  01 

-0.600810-02 

J 

X 

A 

K 

X 

5 

A* 

0.189900-OA 

0.  800020-0 A 

J 

X 

5 

K 

X 

l 

A* 

-0.28552D-02 

0. AA7270-0A 

J 

X 

5 

K 

= 

2 

A* 

-0. A5A27D-03 

-0.809120-02 

J 

X 

5 

K 

X 

3 

A* 

0.565A7D-02 

0.922370-03 

J 

* 

5 

K 

* 

A 

A* 

0.120200-02 

0. 12A380-01 

J 

= 

5 

K, 

X 

5 

A* 

0.997360  00 

-0.  172  7oi3-02 

A 

= 

TOTAL 

CONVERSION 

COEFH  C I ENTS  SLAB  NUMBER  = 9 

J 

X 

l 

K 

X 

1 

A* 

0.130000  01 

0.0 

J 

X 

l 

K 

z 

2 

A* 

0.0 

0.0 

J 

X 

1 

K 

X 

3 

A* 

0.0 

0.0 

J 

X 

1 

K 

X 

A 

A* 

0.0 

0.0 

L 


J = 1 K * 5 
J = 2 K * 1 
J * 2 K = 2 
J = 2 K = 3 
J « 2 K = 4 
J » 2 1C  * 5 
J = 3 K.  = 1 
J * 3 K = 2 
J * 3 K * 3 
J * 3 K = 4 
J - 3 K = 5 
J * 4 K = 1 
J = 4 K = 2 
J = 4 K = 3 
J = 4 K = 4 
J = 4 K = 5 
J = 5 K = 1 
J * b K * 2 
J * 5 K = 3 
J * 5 K ~ 4 
J * 5 K = 5 


I 

I 


A*  0.3  0.0 

A*  0.0  0.0 

A*  0.100000  01  0.0 

As  0.0  0.0 

A*  0.0  0.0 

A*  0.0  0.0 

A = 0.0  0.0 

A*  0.0  0.0 

A=  0.100000  01  0.0 

A=  0.0  0.0 

A*  0.0  0.0 

A*  0.0  0.0 

A*  0.0  0.0 

A*  0.0  0.0 

As  3.100000  01  0.0 

As  0.0  0.0 

A=  0.0  0.0 

A=  0.0  0.0 

As  0.0  0.0 

A*  0.0  0.0 

As  0.100000  01  0.0 
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ELECTRIC  FIELO  STRENGTH  AS  A 

GAMMAI OEGi»  0.0  PHHOEGI* 

l 

•1MOUM)  AMP  1 00 1 ANGtGEGI 

FUNCTION  OF  RHO 

0.0  2TUNI> 

RHGIKM1 

5.000  2RIKN1«  10.000 

X 

AMP  (08)  ANGIOEGI 

RHOIKMI 

r 

AMP ( DB ) 

ANGIOEG) 

25.00 

76.29234 

263.2688 

25.03 

39.90491 

152.2*87 

*5.00 

238.2237 

SC. 00 

72.98376 

257.5801 

50.03 

56.25769 

135. 7691 

50.00 

37.85649 

212.8739 

7S.00 

70.78876 

252.3931 

75.00 

53. 75733 

119.4766 

75. OC 

35.  51668 

187.0150 

1 00.00 

68.98273 

247.8470 

100.00 

51.6l7b9 

103. 3704 

100.00 

33.77069 

161.0755 

125.00 

67.35411 

244.0912 

125.00 

49. 

8 7.  3927 

125.00 

32.39372 

135.1018 

1 S0.00 

65.82755 

241. *690 

150.00 

47.32316 

71.46 8* 

150.00 

31.28989 

109.1563 

1 7b.  00 

64.38651 

239.4891 

175.00 

43.32722 

55.4832 

l 75.00 

30.41103 

83.3227 

2 00.00 

63.04  761 

238.  7 764 

200.00 

42.8  8458 

39.2385 

200.00 

29.72945 

57.7085 

225.00 

61.84  721 

239.0170 

11'*.  00 

40.02911 

*3.02*1 

225.00 

29. 16455 

33.4116 

250.00 

60.81107 

239.8885 

*50.00 

36.53989 

66  Ot> 

250.00 

28.83185 

8.5596 

275.00 

59.95270 

240.9785 

275.00 

31.84071 

341.54  00 

275.00 

28.64261 

344.3625 

TOO. 00 

59.25110 

241.8292 

300.  00 

2 5.  12  791 

297.  7227 

300.00 

28.57390 

320.9485 

325.00 

58.65680 

242.0927 

32  5.00 

24.67183 

209.2661 

325.00 

28.60027 

298.4187 

3 50.00 

58. 10954 

241.5  722 

350.00 

29.65506 

loO»042* 

350.00 

28.57823 

276.8580 

375.00 

57.54008 

240.5263 

375.00 

32.82829 

144.7204 

375.00 

28.73222 

256.3191 

400.00 

56.90408 

238.4413 

400.00 

34.  76332 

127.9352 

400.00 

28.90077 

236.7713 

425.00 

56. 16014 

236.0790 

425.00 

35.99524 

1 i2 .9103 

425.00 

29.06299 

218.1858 

450.00 

55.27319 

233.4316 

450.03 

36.  76811 

98. 8009 

450.00 

29.20108 

200.5223 

475.00 

54.20961 

230.6798 

475. GO 

37.23258 

84.8595 

475.00 

29.28441 

183.4644 

500.00 

52.9315b 

228.1687 

500. 00 

37.36682 

71.8263 

500.00 

29. 33116 

167.5499 

525.00 

51.39867 

226.2026 

52  5.  00 

37.31810 

59. 1317 

525.00 

29.31  720 

152.4182 

550.00 

49.55347 

225.3099 

550.00 

37.08209 

46.7381 

550.00 

29.23508 

138.0355 

575.00 

47.32220 

2 2 6 . 4329 

575.00 

36.67429 

34.6322 

575.00 

29.07872 

124.3776 

600.00 

44.63403 

231.4045 

600.00 

36.13287 

22. 7550 

600.00 

28.88341 

111.2962 

625.00 

41.63359 

244.3454 

62  5.  00 

35.397&1 

11.2109 

625.00 

28.5588  2 

99.0160 

650.00 

39.46725 

269.7253 

650.00 

34.4989b 

0.0290 

650.00 

28.14732 

87.4632 

675.00 

39.95416 

298.2  769 

675.00 

33.42958 

349.2854 

675.00 

27.64709 

76.6715 

700.00 

41.94011 

315.4939 

700.00 

32.17734 

339. 1038 

700.00 

27.05814 

66.6955 

725.00 

43.81900 

323.6453 

725.00 

30.  79295 

330.0586 

725.00 

26.44009 

57.7078 

750.00 

45.29208 

327.J501 

750.00 

29.  14931 

3*1.5710 

750.00 

25.68810 

49.5037 

775.00 

46.39511 

329.0024 

775.00 

2 7.26694 

314.6301 

775.00 

24.86  755 

42.4046 

800.00 

47.21843 

329.7319 

800.00 

25.14587 

310. 1309 

800.00 

24.00182 

36.5465 

825.00 

47.83922 

330.0955 

82  5.00 

22.87070 

309. 5* 7o 

825.00 

23. 12830 

32.0445 

8 50  . 00 

48.31659 

330.4138 

850.00 

*1.31863 

315. 1108 

850.00 

22.36522 

29.1504 

875.00 

48.70630 

330.7258 

875.03 

19.80577 

323.9875 

875.00 

21.65184 

27.1271 

900.00 

49.03757 

331.0835 

900.00 

19.75130 

332.7961 

900.00 

*1.09418 

26.0126 

925.00 

49.33565 

331.4578 

92  5.00 

20.46809 

336.9707 

925.00 

20. 72151 

25.3085 

950.00 

49.61635 

331.7844 

950.00 

21.36862 

336.0334 

950.00 

20.52293 

24.4713 

975.00 

49.88675 

331.9897 

975.00 

*2.156  94 

331.3999 

975.00 

20.45189 

23.0804 

1000.00 

50.  14691 

332.0081 

1000.00 

22.75043 

324. 3b45 

1000.00 

20.44556 

20.9269 

1025.00 

50.39211 

331.7952 

102  5.00 

*3.14766 

315. 7686 

1025.00 

*0.44417 

17.9951 

1050.00 

50.61794 

331.3523 

1050.00 

23.42491 

306.  0242 

1050.00 

20.44395 

14.5671 

1075.00 

50.81192 

330.641 1 

1075.00 

23.48744 

295.5928 

1075.00 

20.33031 

10.3539 

1100.00 

50.96849 

329.7109 

1 100.00 

23.42499 

284.0973 

1103.00 

20. 1 1635 

5.8188 

1125.00 

51.08234 

328.6028 

1125.00 

23.25829 

273.4995 

1125.00 

19. 78918 

1.1468 

1150.00 

51.  15012 

327.3650 

1 150.00 

23.00464 

26*. 1365 

1150.00 

19. 34108 

356.5193 

1175.00 

51.17041 

326.0486 

1175.00 

22.67813 

250.7361 

1175.00 

18. 76810 

352.1216 

1200.00 

51.14354 

324.7034 

1200.00 

22.28987 

239.42  54 

1200.30 

18.07005 

346.1531 

1225.00 

51.07132 

323.3  752 

122  5.00 

21.84828 

228.3357 

1225.00 

17.25197 

344.8381 

| 

1250.00 

50.95685 

322.1050 

1250.03 

21.3  5965 

217.6042 

1250.00 

16.32759 

342.4343 

1275.00 

50.80411 

320.9265 

1275.00 

20.82889 

207.3739 

12  75.00 

15.  32566 

341.2273 

1500.00 

50.61781 

319.8655 

1300.00 

*0.2  610* 

197.7929 

1303.00 

14.29841 

341.4915 

1 

1325.00 

50.40297 

318.9390 

1325.00 

19.66304 

189.0123 

1325.00 

13.32873 

343.3853 

1350.00 

50.16458 

318.1558 

1350.00 

19.04646 

181.1795 

1350.00 

12.52348 

346.7732 

| 
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11  71.00 

49.90  724 

317.1119 

1375.00 

18.43004 

174.4257 

1375.00 

11.97878 

351.0781 

1400.00 

49.63486 

317.3127 

1400.00 

1 7.84206 

Lo8 • 8399 

1400.00 

11. 72717 

355.3843 

1425.30 

49.31028 

316.6341 

1425.00 

17.  3202  5 

164.42 94 

1425.00 

11.71518 

358.8298 

1410.00 

49.01119 

316.3669 

1450.00 

16.90714 

lot.0745 

1450.00 

11.63764 

3.9493 

1471.00 

48. 71000 

316.1960 

1475.00 

16.63916 

138.5040 

1475.00 

11.96946 

1.7006 

1100.00 

48.43399 

316.1116 

1500.00 

16.53325 

130.3272 

1500.00 

12.09305 

1.2926 

1121.00 

48. 10119 

316.1091 

1525.00 

16.57874 

1J4.12  79 

1525.00 

12.10005 

3.0219 

1110.00 

47. 76273 

316.1919 

1550.00 

16. 74065 

131.371b 

1550.00 

11.98337 

356.1875 

11  71.00 

47.40  327 

316.3713 

1575.00 

16.97197 

ltd . 4663 

1575.00 

11.72  742 

356.0647 

loOO.OO 

47.02144 

316.6677 

1600.00 

17.22643 

I4t. 7389 

1600.00 
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39.9071  7 

07.  1379 

3575.00 

-6.90990 

192.0618 

3575.00 

0.59286 

130.6951 

5o  00. 00 

39. 88676 

bd.2501 

3600. 00 

-7.91166 

199.  3673 

3600.00 

0. 5o586 

130.1632 

3o25.00 

39.66  390 

65.3131 

3625.00 

-8. 9o**  73 

197.  1*06 

3625.00 

0.53613 

129. 7257 

30  50. 0 0 

39.83633 

70. 5313 

3c50. 00 

-10.09030 

lal.5008 

3650.00 

0.50361 

129.3257 

3675.00 

39.80998 

71.5096 

3675.00 

-11.08370 

15  7. 9906 

3675.00 

0.9ud89 

12o  .9570 

3700.00 

39.76881 

72. *530 

3 700.00 

-12.00968 

lu5.2bl 8 

3700.00 

0.96115 

128.6138 

3725.00 

39.72796 

73.1 669 

3 72  5.00 

— 12.6  82 1 6 

1 79. 6o02 

3725.00 

9.39029 

1*6.2913 

3750.00 

39.68063 

79.0551 

3 750.00 

-13.00503 

165.0171 

3750.00 

0.39569 

127.9859 

3775.00 

39.6*819 

79.5290 

3775.00 

-12.93582 

195.2635 

3775.00 

0. 2 96  To 

127.6951 

3800.00 

39.5698  7 

75.  7 779 

3600. 03 

-12.53720 

2 39.93  69 

3800.00 

0.2*291 

127.9177 

3825.00 

39.50578 

76.6215 

3825.03 

-11.92770 

211.9312 

3825.CC 

0. 163.1 

127.1539 

3850.00 

39.93  591 

77.5596 

3850.  00 

-11.22000 

21  7.  7959 

3850.00 

0.11 7o0 

126.9056 

38  75.00 

39.36039 

76.2958 

3&75.00 

-10.99221 

222.2127 

3875.00 

0.09969 

120.6757 

3900.00 

39.27922 

79.1350 

3900.00 

-9.78909 

2*5. 962  7 

3900.00 

-0.03515 

126.9689 

3925.00 

39.19275 

79.9d0  7 

3925.00 

-9.  132  79 

2*7.6016 

3925.00 

-0. 12296 

1*6.2892 

3950.00 

39.10117 

60.0368 

3950.00 

- d. a32 1 6 

229.9371 

3950. OC 

-0.21353 

126.199'' 

39  75.00 

39.009  76 

81.7067 

3975. 00 

-7.98961 

230.3280 

3975.00 

-0.32182 

120.0399 

9000.00 

33.90387 

82.5936 

9000.00 

-7.50390 

261.  1999 

9000.00 

-0.93251 

125.9823 

9025.00 

33.79887 

83.5005 

9C25.00 

-7.37290 

231.3293 

9025.00 

-C. 55006 

125.9793 

9050.00 

33.69019 

89.9301 

9050. 00 

-6.69229 

26*  *o0 1* 

*050.00 

-0.67367 

126.0305 

90  75.00 

33.57826 

85.3699 

9075.00 

-6.  36012 

2 3 1 • 9o  5 1 

9075.00 

-0.80235 

12o • l 595 

9100.00 

33.96359 

86.3668 

9100.03 

-6.37313 

231.1697 

9100.00 

-0.93990 

12b. 3525 

9125.00 

33.39668 

8 7.3  77o 

9125.00 

-5.82850 

2a0. 7399 

9125.00 

-1.07030 

126.&189 

91  50.00 

33.22610 

88.9166 

9150.00 

-5.62390 

230.2050 

9150.00 

-1.20618 

1 2o .9580 

91  75.00 

33.10890 

89.6907 

9175.00 

- 5.  95o8  7 

2 * 9 . 39  9 3 

91 75.00 

-1. 39199 

127.3727 

9200.00 

32.98819 

90.5953 

9*00.00 

-5.32539 

2*8. 9385 

9200.00 

-1.97576 

127.8589 

9225.00 

32.86  792 

91.7293 

9225.30 

-5.22752 

2*8.2905 

9225.00 

-1.60620 

128.9117 

92  50.00 

32. 79  82  9 

92.6953 

9*50.00 

-5.  16151 

2*7.5*07 

*250.00 

-1.73192 

129.0267 

92  75.00 

32.62981 

95.0915 

9275.00 

-5.  12572 

226. 7927 

92  75.00 

-1.85176 

129 .6962 

9300.00 

32.51302 

95.5163 

9300.00 

-5.11869 

2*0. 0689 

9300.00 

-1.96979 

130.9117 

9325.00 

32.39891 

96.5679 

932  5.00 

-5.10890 

2*5.3583 

9325.00 

-2.07031 

131 . 1637 

9350.00 

32.28692 

97.6663 

9350.00 

-5.18526 

*29.0719 

9350.00 

-2.  16791 

131.9929 

43  75. 00 

32.17747 

99.1429 

4375.00 

-5.25658 

224.0178 

4375.00 

-2.25746 

132.7377 

54  00.00 

32.07191 

1.00.4608 

4400.00 

-5.35164 

223.4043 

4400.00 

-2.33907 

133.5397 

5425.00 

31.97002 

101.7950 

4425.00 

-5.47015 

222.8391 

4425.03 

-2.41311 

134.3393 

44  50.00 

31.87202 

103.1423 

4450.00 

-5.61070 

222.3299 

4450.00 

-2.48018 

135. 1283 

44  75.00 

31.77806 

104.4992 

4475.00 

-5.77276 

221.8845 

44  75.00 

-2.54101 

135.8992 

4500.00 

31.6882  3 

105.8623 

4500.  00 

-5.95564 

221. 51 12 

4500.00 

-2.59651 

136.6463 

4525.00 

31.60255 

107.2283 

4525.00 

-6.15870 

221.2188 

4525.00 

-2.64  765 

137.3646 

4550.00 

31.52095 

108.5938 

4550.00 

-0.38128 

221.0166 

4550.00 

-2.69544 

138.0508 

45  75.00 

31.44331 

109.9555 

4575.  00 

-6.  62266 

223.9152 

45  75.00 

-2. 74091 

138.7028 

4600.00 

31.36946 

111.3105 

4600.  00 

-6.862  00 

220.9255 

4600.00 

-2.78508 

139.3194 

4625.00 

31.29918 

112.6559 

4625.00 

-7.  15834 

221.0598 

4625.00 

-2.82891 

139.9006 

4650.00 

31.23216 

113.9893 

4650.00 

-7.45047 

221.3310 

4650.00 

-2.87328 

140.4474 

46  75.00 

31.16811 

115.3086 

4675.00 

-7.  75691 

221.7527 

46  75.00 

-2.91901 

140.9613 

4700.00 

31. 10664 

116.6118 

4 700.  00 

-8.07582 

222.3390 

4700.00 

-2.96681 

141.4446 

4725.00 

31.04741 

117.6975 

4 725.  00 

-8.40497 

223.  1041 

4725.00 

-3.01 730 

141.8994 

4750.00 

30.99002 

119.1645 

4750.00 

-8. 74161 

224.0617 

4750.00 

-3.07099 

142.3304 

47  75.00 

30.93405 

120.4122 

4775.00 

-9.0  62  49 

225.2244 

4775.00 

-3.  12830 

142. 7392 

4600.00 

30.87914 

121.6402 

4800.00 

-9.42377 

226.6025 

4800.00 

-3. 18956 

143.1295 

4825.00 

30.82486 

122.b483 

4825.00 

-9. 76103 

228.2032 

4825.00 

-3.25500 

143.5047 

4850.00 

30. 77087 

124.0368 

4850. 00 

-13.03935 

230.02  94 

4850.00 

-3.32479 

145.8679 

48  75.00 

30.71670 

125.2063 

4675.00 

-10.4333o 

232.0780 

4875.00 

-3. 39903 

144.2221 

4900.00 

30.66229 

1 26.3  57o 

4900.00 

-10.69747 

254.3387 

4900.00 

-3.47774 

144.5701 

4925.00 

30.60  707 

127.491  7 

492  5.00 

- 10. 966 1 5 

236.  7935 

4925.00 

-3.56092 

144.9146 

4950.00 

30.55086 

128.0099 

4950.00 

-11.20425 

239.4155 

4950.00 

-3.  64855 

145.2577 

49  75.00 

30.49339 

129.7134 

4975.  00 

-1  1.40738 

242.1700 

4975.00 

-3.  74054 

145.6017 

5000.00 

30.4344a 

130.6039 

5000. 00 

-11.5  7233 

245.0155 

5000.00 

-3. 83684 

145.9482 

IV.  Program  Checks  and  Some  Results 

The  sample  input  output  case  discussed  in  the  previous  section  is  an 
example  of  one  of  a number  of  checks  made  of  the  current  program.  That  par- 
ticular case  represents  propagation  in  a very  slightly  inhomogeneous  guide  in 
which  the  horizontal  inhomogeneity  begins  just  25  km  from  the  transmitter  (i.e. 
XVAL(8)  = 25.0  km  since  a nine  slab  model  was  used  for  the  inhomogeneity)  and 
is  1000  km  in  extent.  The  ionosphere  is  described  by  exponentials  and  in  the 
convention  of  Wait  and  Spies^varies  linearly  over  its  1000  km  extent  from  H1  = 

70  km  to  H'  = 71  km  with  constant  0 = 0.3  km’1.  This  is  a sufficiently  modest 
inhomogeneity  that  one  would  expect  the  mode  conversion  and  WKB  results  to  be 
identical  for  all  practical  purposes.  Figures  3 through  5 show  the  mode  con- 
version and  WKB  results  for  the  electric  field  components  E , E and  E for 

z x y 

four  orientations  of  the  transmitter.  The  latter  was  at  5 km  and  the  receiver 
at  10  km.  Examination  of  the  plots  will  show  that  the  mode  conversion  and  WKB 
results  are  indeed  very  nearly  identical. 

As  an  example  of  a more  realistic  terminator  problem.  Figures  6 through  8 
show  the  Ez,  Ex  and  Ey  results  for  a Hawaii-San  Diego  path  as  a function  of 
terminator  location  (i.e.  the  abscissa  is  the  distance  from  the  transmitter  to 
XVAL  (20)  since  21  slabs  were  used  to  model  the  terminator)  for  a transmitter  and 
receiver  altitude  of  15  km.  The  terminator  extends  1000  km  and  the  ionosphere 
(again  assumed  exponential)  varies  from  H1  = 86  km,  & = 0.5  km’1  to  H'  = 70  km, 

6 = 0.3  km’1.  Curves  for  four  orientations  of  the  transmitter  are  shown  on 
each  plot.  Some  jaggedness  in  some  of  the  plots  occurs  as  new  slabs  pass  over 
the  transmitter  or  receiver.  However,  the  jaggedness  would  have  been  much  worse 
for  the  Ey  case  except  for  the  fact  that  the  fundamental  input  from  the  wave- 
guide runs  (i.e.  the  program  of  reference  4 ) was  generated  using  tolerances 

-5 

of  10  degrees  for  both  real  and  imaginary  partsof  the  eigenangle  as  compared 
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with  our  usual  tolerances  of  0.01°  for  the  real  part  and  0.005°  for  the 
imaginary  part.  The  increased  tolerance  results  in  only  a slight  increase  in 
execution  time  of  the  waveguide  program  and  we  would  recommend  using  the  smaller 
tolerance  when  generating  data  for  the  present  program  and  when  fields  are 
of  interest. 
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APPENDIX: 

LISTING  OF  PROGRAM 

C SIMPLIFIED  MODE  CONVERSION  MODEL  MODIFIED  TJ  CALCULATE  PIcLDS 
C PUK  AN  ANTENNa  OP  ARBITkAKY  HEIGHT  AND  ORIENTATION. 

IMPLICIT  REAL  *81A-H,U-Z) 

CUMMON/HGT  EMP/FP1I25*  'j)  * PP  A(  23 » :> ) 

COMMON/ 1 ERM/WT ,NTR 

COMMI  N/CAP/CAPI  125,5,5)  . TNQRMI  3) 

CUM MON/NL  INPT/ThtTA (<.:>,  5) , P'LFR  I 2 a , 5 ) , XT  KA  l 3 , 3 , 2 j,  5 ) ,TUPHT  (2  5)  , 
i X V A L lit,)  ,FPE'.',KHOMAX»RHOMlN,DELRhO,DELTAX,EPSLON(2:>)  , 

% SIGMAI2.,)  ,NRSLAB,NRMUDE  ,NTMAX 

IdMM  N/MCSTOk/  A 1 2 t» , , 5)  , SI  2:>,  5)  ,C  (25,  5)  ,NTHSG(25)  ,KVR  AuT.nvRATT  , 

$ AVRK.OT  iAVRKTT  , CONST, OMEGA, WAVE  NO 

DUMMbN/MCFLUT/FUOO  ) , Di3l  J , , 400)  , ANG  I i,  H»  <*Q  J ) , I UPLUT  ( 131,1  SUb.NRP 
GUPiMON/X PL  C T/  AM  IN  , X I N C , V Ml  N,  Y I NC  , S I ZEX  , SI  ZE  Y 
L UMML  N/Ho  I nPT/G  A MM  A ( H } ,Phi(4)»ZT  ,ZR 
t , SING  AMI  A)  ,CUSGAM(4J  ,SINPHI  141  ,CL)SPHl(41 
DIMENSION  EUFPE R ( 2u00 ) 

DIMENSION  Z ( 2 1 

HEAt*<«  bCDl  20  ) , Out/ • LE  N ' / 

K F A L * 4 R,DL',AI.G,XMIN,XINC,YMIN,Y1NC  ,SUEX,SIZEY 
CuMPLtX*lt  PP1,PP3 

COM  PI FX*10  THETA, A, S, C, PUP R,  IM/ (O.OLO, I . OUO I / , C AP I , TN ORM 
LUMP!  EX*  it  XTRA,NTHSw,T  l ,T2,  T3 » T* 

DATA  TkOP 1/6. 2biIB5D0/, VEL I TE/2. 49792805/ , ALPHA /3. I4D-4/ , 

& UEokAD/ l .7453240-^/ 
uAT  A LAST/O/,  IPRNTA/C/ 
kF  AL*o  R VP  ACT , K VR  AT  T 

named i st/ da  tom/  rhomax, 

L RmUMIN.UElT A A, NR SL A B , NK MODE , NT  MAX » XV AL , r Rt w » I OPLuT , t P SL jN, 

1 SIGMA.C  E LRHC, I P I RS T , L AS T , I PL T OP , 

$ XM  IN , X I NC ,YMIN,Y INC, SIZEX  ,SIZEY , 

$ GAMMA, PHI, ZT,ZR,NRP, IPRNTA 
i , I NTP  LG 
C 

1 NTP  lg  = o 

IuPLOT(l)  = 0 
10  DO  250  K= 1 , 2 00 
REAL! 5,201)  BCD 
*P I T E I o , 2G2 ) BCD 
rip  1 TE( i , 2C1 1 BCD 
IFlbCD(I)  .EQ.  OOT)  GO  TO  260 
250  CONTINOE 
260  REWIND  I 

‘•EACI  1, DATUM) 

REWIND  1 

00  200  N=I,NRP 

oAMMAIN)  * GAMMA(N)*DEGRAD 

PH]  IN)  = PHI  IN) *OEGP  AD 

S:\oAMlN)  = OSINIGAMMAIN) ) 

CuSoAM(N)  = JCOSIGAMMA(N) ) 

SINPhilN)  = OS  INI  PH  1 1 N)  ) 

CLSPhl IN)  = DCOSIPH  l(N)  ) 

200  CONTINUE 
PRINT  i 02 


52 


I J? 

iu- 
l 0^ 
iP- 


1 ID 


1 3D 


i 35 


DO  11'  M-l  ,NRSLAj 

i > i iU  K - 1 , M-  MiJt'E 

'■rAt,  1 uo  f TFirTA{M»K)»Ti,T2 

si;Ai.  100,  T HE  7 A ( M t K \ * T3  » Ta 

• • V ? 1 , r r'F»  ( »7LlPriT  (.-IJ 

^1  mI  ±03,  TFlET  AIM  ,K  ) , T l ,T2,  T3,  ,FuFk(M,iO  , T GPHT  ( MJ 
I V.RS,.  it  5A,  • THtTA*  ,13A  ,*  T l*  ,2Ja,  • • ,20X,*  T3*  , 20  X , ' U*  , . Da  , ' Fuf  r • , 

IPX,  • fOPhT*  ) 

rtK/i'  T ( 1 * »i’F7.it4.X»  5(4Di.O»i,4X),F<,,li 

I UK  MAI  ( lA»2F'*.3,iX,‘,ui3.o) 


F c k *■  A T | •tiXttiUo  • o , 4.  X , F !)«  «.  1 
')(v,nJ  = LDSIMTHf  TAiM,M*Ot&KAD) 
t„ I M , i\i  = L DC US  I Tri  t T A ( K • K j *uE bk  AD  ) 


XT3  A(  1 , i ,y  , K) 
aTR A I l,2,M,XJ 
X TF'  A { 1 t 3 i M i K) 
A T R A ( , 1 » M * K ) 

aTRAI2,2,K,K> 
J T h a ( i j t K t K ) 
,» T A A ( 3 , 4,1-tM 
aT*  A(  3 , 4 , M , K ) 
XTR  M 3,3,K  ,K) 


= T1*SIM,a)**2 
= Tl*SlM,K) 

= -T3*3(M,M 

= - T 1 
= 1 3 

= — T j *•  7 t*  S I f : t ft ) 
= -T3*I  <t 
= T 2 


CLM  I ME 

iKIlKCl  ID  .UF  . 0 • A,mD.  IF  IKS  I 


i\  E . 


OCALL  PLOTSI  bUFi-ck, 4030,1  j) 


N T — 1 

*AVE  Nj  = T*U  FI  *100u.J*FkEw/VEL  IT  E 
v-UNSl  = C.C324d*hAvt  M 0/ OS  CRT  I F R E w ) 
liME  oa  — 7 iaL  P 1 * F KE  0 *1  DO  0 • 
ft  V * ;i  C T = LL/P(  DLJbl  *aVT  NO/  ALPfi  A J / 3 . ) 
A»ki,Il  - RvRAOT#*2 
a V R ft  c T = 1 . / K V K Ai ) T 

A '» R ft  J T = A VRKOT  **2*0.  P 
Do  i j 0 L=  1 « UR  SL  Ab 
uTHSCIL)  = L .+ALPHA*1 OPHT l Li 
ll-Ll-  = 0 

00  l^b  M-  1 , NR  SL  Ab 
v.K  = ,\KSIA|,-M+1 

1 F ( MM  . Nt  • it|R  Sl  A J ) IF  Lb  = 1 

CULL  HT  IImTL  (CAP  1 , TIMOR  M,  lFLb,  MM,  I NT  FIG) 


All)  = ZT 
Z i 4 ) = ZK 
CALL  FT(jA  IN  ( Z 1 
00  1 3 o M=  1 , FjR  SL  A rt 
00  130  K=1,ImWMO0E 
ATRA(ifltMfK)  = X T K A l 1 , 
XTR  A(  1 , 2 »M  »K ) = XTRAl  1, 
XTR  AIl»3,M,Kl  = X T K A ( 1 , 
\TRi*li_,l,M,K)  = XTRA<  4 , 
XTPA(2,2,M,k)  = X TR  A l Z • 

XTRA(Z,3,M,K)  = XTR  At.  Z, 
XTRAIj,I,M,K)  = X T P A t 3 , 
XTR  X. ( 3,t,M»k)  = X1PA(3, 

Alft  A| 3, 3,M,Ki  - XTK Al 3, 


1 , M , K I * FF  l ( M 

2 ,M,M*FFl  l M 

3 , M , ft  I * F F 1 I M 
l ,M,K)*FF  1 IM 
4,M,K»MFFilM 
3,M,a)*FF  1 IM 
1,M,K)*FF1(M 

4 ,M  , K ) *FF  i ( Vi 
3 * M , K ) * F r 3 i M 


» K ) **Z 
,K )**2 

, M * F F 3 i M , K)  /FUFk(  M,M 

,K i**Z 

,k)**2 

,k>*FF3i  M,M  /F0FR(  M,XJ 
,K)*FF3(M,k)/FOFk(M,Kl 
,K)*FF3{M^)/FL:Fk(M,M 
,n)**Z  / (F0FR(M,KJ  + 't<-J 
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1 3t»  CGN  T 1 NOE 

113  IF  ( XVAL(NRSLAb-i)  .Gt.  0.1  GO  TO  ill 


'Jit 

112 

J= 1 » NRMuDE 

00 

112 

K=l , NRMUJF 

Jl; 

1U 

L = 1 ,NRSL AB 

« ! L 

t K t 

JJ  = 0.0 

DC 

113 

L-2 .NRSLAB 

IF(XVAL(NK$LAE*l-L)  .GE.  O.J  GO  TO  il9 
113  CONTINUE 
NT*  = 1 
GO  TO  il  / 

lit  NTK  = NKSLAB+2-L 
1 1 7 CONTINUE 

00  ilo  J= 1 i NRMODE 
DO  116  K= 1 , NR  MODE 

lib  IE(n  .bO.  J)  A!  NTR  , K,  J)  = ( 1. 0.0.0) 

NR MO  = NRSLAB-1 
JO  4C1  MM* l, NTR 
M = NTR-MM+l 
♦ 01  CALL  MUST  EP l M ) 

1 E ( 1PKNTA  ,EU.  0)  GO  TO  9 1 
PRINT  905 
DO  451  L=1 , NTR 
PRINT  900, L 
DC  t5 1 J=  1 , NRMOOE 
DO  451  K = 1 » NRMODE 
PRINT  901.  J,K,A(L,  J.N) 
t5i  CONTINUE 

91  IF( IPLTOP  ,E0.  U CAlL  MCFLD 
IF! IPLTUP  .Ew.  2)  CALL  MCFLD2 
NT  = NT* l 
JU  1C  o ME= 1 . NRMU 
10o  XVAL(ME)  = XVAL ( ME  J ♦D  EL  T AX 

IFIXVAL(NTR)  .GE.  0.  .AND.  NT  .LE.  NT  MAX ) Go  TO  118 
1FINT  .LE.  NTMAX)  oU  TO  91 
l F (LAST  . L«.  0)  GU  TO  10 

I F ( 1 0 PL  OT ( 1 ) .NE.  0)  CALL  PL OT ( 0 . , 0 . , 999 J 
HETUPN 

til  NTK  = NRSLAB 

DO  150  J = l.  NRMODE 
DO  i.  50  K * 1»  NRMODE 
DO  120  L = 1, NRSLAB 
A ( L , K , J ) = 0.0 
120  CONTINUE 

IFIK  .EU.  J)  AINRSLAb.K. J)  = (l. 0,0.0) 

1 30  CONTINUE 

NR  MO  = NRSLAB  - 1 
c THE  LOOP  400  DETERMINES! A) 

C IN  SUCCESSIVE  SLAbS. 

JO  4 00  MM  = 1,  NRSLAB 
M = NRSLAB  - MM  * 1 
CALL  MCSTEP(M) 

IFU.xSLAb  .LE.  1)  RETURN 
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40J  Cl.VJINUE 

II  l ll-RNTA  ,'rj.  0)  Gu  T 0 90 

J'MM  iQj 

Of.  9 SO  L = 1 t NR  SLA? 

P K 1 i\  1 vOOtl 
u0  -,bJ  J= 1 , NRMuJE 
Mj  <til  K=l*MMOJc 
P l»  1 N T 9 0i,J«K,A(L,J,lv) 

'♦b  J CONTINUE 

9f  IK  IPLTOF  .EU.  1)  call  i-lCr-LU 

i I-  ( I H L T UP  .bv.  2)  CAlL  KCFlD2 
NT  = NT  ♦ l 

JO  iOi)  ML-  = 1»  ImRMO 

A V A L ( ME ) = XVAL(ME)  ■*•  OELTAX 
103  CONTINUE 

IMNT  .Lt.  MMAXl  GO  TO  90 
IFILAST  • t (j  . 0)  Gu  Tl  10 

l F ( 1 D PLOT  ( 1 ) ..ME.  0)  CALL  PLOT (C ., 0 ., 999) 

•'Ll  UP  .v 

201  R>*ATl2^A4) 

*0*  hf  «MaT{  • • ,20 A*,) 

900  EuPMATIlh  , 1*X, 

i *A  = TOTAL  CONVERSION  COEFFICIENTS • » oX  » * SL  Afl  NUMtiEk  = *,lc,/) 
9Ci  FOMlU14A,'  j =',d,3X,'  « = • * !2»bX,'  A* • , ( c 1 b.  b » cl  b . b ) , / ) 
90b  FORM  AT  1 1H1 ) 
t NO 


) 
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SUBklJUT  INC  HTINTLICAp I .NORM,  IFLG.M,  INTFLG) 

C CALCULATE  normal  I/AT  ION  INTEGRALS  AND  INTEGRALS  OF  riEIonT  j A 1 N S IN 
L AJJmCENT  iLABS. 

IMPLICIT  REAL  *3(A-H,U-Z) 

CuMrtLN/MCiM‘T/TMETAl^:>»a)»FOFK(2:>».>)»xTRAt.>»3»2.>»:>)»TJPhT(2:>)» 

> XV  ALUS)  ,MEw,RHU'1Ax,KHuMlN,UtLRhO»JtLTAX,EPSL'JN(2.>J  , 

% SIGMA(2!>)  .NRSLAB,  NKMODE  ,NTMAX 

lI'KMlN/MCS  I OR/ A (2  3,  t> , 5)  ,S(2S,p),Ct25,i>l  ,NThSuI2  3i  ,nVR  AJT  , nVK  ATT  , 

* AVKXQT ,AVKKTT .CONST, OMEGA, aAVfc NO 

CGMfLtX*lo  NTHSC 

oOMPL  tX*ib  PThA.HlTA , H2  TA  .HlPRTA  ,H2PKTA  , Hr  THAI  5 ) , EYTrtAI  Dj  , 

$ HYTHPAIB) , eythpai t>  j 

CIjMPLEX*U  THETA,  FoFR,  a,  S,  C,  SSJ,  CSU,  IM/<  0.  JO,  l.DO) /,WGSm, 

L SwKUjT,kTIOKT,PO.PTH,HlO,H2O,HlPRM0,H2PRMO,CAPHlJ, CAPH20, 

{.  A1ST,«2NC,A3RD,A4Th,0£n12,  DEN  34,  DENMF  ,NUKMF  , 

L H1T.H2T ,HlPRMT,H2PkMT ,HYTH(3) , EY Th ( a ) , HY THPk ( 3 i , E Y TnPR(b), 

L HYOPKI  b)  , E YOpR  Cb)»EYD(5)  , MULT  , F AC  l , F AC2,  NORM  I 2 a,  3 , S) ,PS(S» , 

$ CAPI(*:>,b,S),PHVlri<5)  , PHYTHP  ( 5 J , PEYTN(  51  » PfcY  TriP  ( ui  » PE Y J(  5 ) , 

L Pf-  YOPR  ( b)  , PH  YOPR  I b ) , XTkA 
kEAL*B  KVRAOT ,K VRATT 
UATA  EPSLN0/6.8b434D- 12/ 

C 

C 

00  100  K = l,  NRMCOF 
iSU  = S ( M , K 1 **2 
CSG  = C I M , iO  **2 

NGSG  = (EPSLONIM)  - lM*SlGMA(rt)/UMEGA J/EPSLNO 
SwkCOT  = LUSORTINGSw  - SSO) 

KSvR  = SwRCiOT 

lFIksak  .LT.  C.)  SQhUOT--SQROOT 

kTIOkT  = L . /NGSU*SOR JOT 

PO  = kVRAtT  *CSU 

PTH  = KVR  ATT*(NThSy. IMl-SSU) 

CALL  MuHNKL IPO, HIO, M2  0, H1PRMJ,H2PRM01 

CAPHiO  = HIPRMO  ♦ AV*KTT*H10 

CAPH20  = h2  PRMO  *■  AVkKTT*H20 

AlSt  = CAPH20  - IM*kT I OR T*KVRAU T*H20 

A2NC  = CAPHIO  - IM*rT10RT*KVRA0T*H10 

A JR  D = H2PFM0  - IM*KVRA0T*SURG0T*H20 

A4TH  = HIPRMO  - I M*NVRAUT*SURJ0T*H1 0 

DEN 12  = H2C*A2ND  - HiO+AlST 

DEN  34  = H2Q*A4TH  - HiO*A:»RO 

CALL  MDHNKl (PTH, HIT ,h2T , HIPKMT ,H2PRMT ) 

HYTh(K.»  = ( H2  T*  A2  NO  - H IT* Al ST ) /DEN  12 
EYTHU)  = (h2T*A4TK  - H 1T*A3RD1 /GENJ4*FUFR < M.Nl 
HYTHPRIlO  = (H2PRMT*A2N0  - H IPRMT *A  IS T J /3EN  12 
tYTHPR(K)  = (H2PRMT*A4TH  - H l PRM T*A3R0 i /DEN34*F OFR ( M , K) 

HYOPklKJ  = I PR  MO*  A^NO  - HI  PRMD*A  l ST  * / DENI  2 
EYOPRIK)  = ItI2PRMC*A<t  TH  - HI  PRMO  *A3PD ) / UEN34*F0r  R I M , K 1 
I F ( IF  LG  .(Y.  C)  GO  TO  100 
P TH  A = KVkATT* l NTHSOt M*1 )-SSO) 

CALL  MDHNKL (PThA.HI T A , H2 T A , H 1 PR T A , H2PRT Al 
HYTHA(K)  = ( H2T  A* A2NU -H 1 T A*A 1ST ) /DEN  12 
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rYThAlK)  = ( H2  T A*  AM  H -H  IT  A*  A 3k  [))  / UE  N3h*  FuFk  l M , R ) 
nVIH  «iK)  = ( FicPRT  A*«4NU-HiPRT  A * A L S T ) / DE  N 1 <L 
lYThFAIKI  = ( H^PRTA*AATn-HiPRTA*A.tkL) /DENjA*FuFK1  M,<0 

1 vJ  J : i J ( n ) = F L Fk l M ,N ) 

iM  Ip'UUi  .fw.  1)  PRINT  l,  Jt.  ,M 
0 ttO  J = 1 1 NKMCJDt 
o.r  2vJ  R = l,MMuDE 
1F(j  .Lu.  F.  I GO  TO  1..0 

•’Cl  - AVRKCT/ l (3 (M, J ) - S(M,M  )**AVENuJ 

t A*.  1 = EY  TH(  R ) * EY  THF  h 1 J ) - E YTH(  J ) ■»  EY  ThPK  { is  ) ♦ nY  Th(n  I^llYTrlPR  to) 

i -HYTHt J i *nY 1 nPR l K) 

FACi  = -EYOlK  )*:YOPP(  J)  ♦ EY  0 1 J I * EY  OP  R ( R I - HYjPk(J)  HYJPMiO 
moi-  iV ( .vi  t J i K i = MUl.  T * ( F At  1 + h AC2  ) 

Ul  IMFLO  .FU.  1)  PRINT  SOS,  M,  J , R , NCRMi  M,  j , is) 

GO  1C.  2hl 

llu  MOL  1 = 2.0*S<H,  JMKVisAUT/rfAYENu 
PTH  = RVkATTMNTHSC'tP  )-SiM,J)**2) 

P"  = RVRATT*C(M,J )**2 

r AC  1 = EYThPK  ( J J **2  t hY  ThPR  < J )**2  ♦ PTH*  ( t yTH(  j)  **2  *■  HYlHlJ)**2) 
F«C  2 = -EYuPRi  J)**c  - FiYOPR  ( - P 0*  l E YO  l J I * *<.  * l.Ji 

N0RM(M,J,K)  - i^UL  T * IF  AC  L*FaC2) 

1F( iNTFLO  .Eo.  LI  PRINT  908 , M,  J , K , A t R M<  M,  j , R) 

2 -to  CON  T 1 NUt 

IF  IIFLG  .to.  0)  GL  TO  SOj 
JC  tOO  K = 1,  NRMODE 
OJ  Ao 0 J = 1,  NXMUDfc 

il-L  I = AvRR JT/ ({ PS( J)  - S(  M » R ) ) * n AV  ElsU ) 
f AC  j.  = EYTHA(K)*PEYThP(  JI-PEYThl  J)*EYTRlPA(is) 

$ FHY  Tn  A ( K J *FFiYThP(  J i-F  HY  TH(  J ) *HYThPA(K  ) 

F AC  t.  = -FVo(KJ*PEYOPF  < J)  * P EY  0 ( J ) * E YORK  I R ) -PUYOPkIJJ  * nY'JPKIK) 
CAP  I ( M,  K , J ) = MULT* (F AC1RFAC2) 

i F C INTEL u .60.  1)  PRINT  9 1 0 , M , R , J , C AP 1 ( M, R , J ) 
m JO  C UN T I NUE 

^00  JO  cJO  J = 1«  NR MODE 
* S l J > = S ( h , J ) 

PHYTR'(J)  = HYTH1  J) 

PhYTFiP(J)  = HY  T HPR  1 Ji 
FEYThtJ)  = EYTH(J) 

PFYlRPt^)  = FYTHPh(J) 

PHYCPPIJ)  = HYOPR(J) 

PEYO(J)  = EYOU) 

PE  Y OPR ( J I = EYOPR(J) 
oOC  CONTINUE 
7cJ  CONTINUE 
RETURN 

9Co  FURMAT (• J* »20X, • INTEGRALS  IN  3LA6',Ia,/) 

JOrt  F0P.MAT(2LX,  'NORM!  ' , Ii  , ' ,'  , II , • , • , 1 1 , • I =*,2013.0) 

9 i J F'UKMAT(21X , *CAPI(  * , II  , * » • , II t * t * , 11 , * ) =',iUij.o) 

END 
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SUBROUTINE  MCFLO 
IMPLICIT  REAL  *o(A-h, o-Z) 

C COMPUTE  f IF. IDS  FROM  AvAL  M I N TO  XVAL  MAX  FJR  T wU  xMTR-kCvk  DISTANCES 
C AT  DLLTAx  INTERVALS 
C 

LIjM  Ml.N  / T Kk  M / i\iT  * N T k 

UMML  N/SPlUT/SAVEU14U2)  , Y 1(  » V2(3,4,402) 

CuMMLN/MC  INPT/ThETAU:>,S),F0FR(2S,b),XTRA(3,.i,2b,b)  , TUPhT  (2  Si  , 

J.  X V AC  ( 2 b)  » FR E w » R HUM  AX  . KHLH 1 N»  Dt cKHO  » DEL  T A X, EP  SL  ON  (2b)  , 

i SIGMA(Zr) ,Nk SLAB, NR MODE, NT  MAX 

CORMLN/MCSTOk/AUb.b,  b)  , S ( 2 b , bi  , C ( 2 S , b)  , NTH Sw ( *. bi  ,nVk  A uT  , xVKATT  , 
i AVkKOT ,A  VRkTT .cons  t , umlga , wa ve  no 

l.OMMl.N/MC  PL  OT/R  (tOB  ) . Oc  l 3 ,4,400)  , AMO ( i , 4,  tO 0 ) , (OPlOT  l 1 J ) , I S Ub , NR P 
CUMHUN/hTGN/M  3 ,2b  ,b,  2) 

CuMM0.'4/HuIMPT/(,A,MMA(4  ) , PHl(4),ZT,ZR 
% ,S INGAM( 4) ,C0SGAM(4i , SINPHI (4) .LUSPHl (4) 

REAL*  4 R.Dt  , ANG,SAVEU,Y1,Y2,  ANG  1(3,4,40  2)  , AnG2(  3,4,402) 

KtAL*d  KVR AOT.KVRATT 

CUMPLfcX*lfc  S0LNA(b,J,4i,IhETA,A,S,C,XTKA,TB,TDBL,TA,F0l-K,r,NTHSw, 

$ IM/U.UDO,  l.JDO)/ 

DATA  Et-AD/o. 37003/ 


(SUB  ^ L 
MP=-i 0 

RHQ  = RHO  KIN 
o'U  CONTINUE 

DU  I LL=2,NKSLAb 

1 F(  XV  AL(NRSLAo«-L-lL  J-RHO  . ob  . 0.)  GO  TO  2 

1 CONTINUE 
M=  1 

Go  T l.  3 

2 M = NRSlAB*2-LL 

3 CONTINUE 

i F ( M . EC.  MP)  Gu  Tl.  72^ 

on  71  0 N=  1 » NRP 

00  71  J L=1,J 

UU  71 C J = l.NPMOOE 

SUL  N A ( J , L , N ) = ( u.G, 0.0) 

no  710  K = L , NRMuDE 

I F ( M .NE,  NTR ) GO  TU  3» 

SOL N A ( J , L , N ) = SULN  A(J,L,N) 

i *A(M,  J,K)*UTRA(  l,L,NTR,K)*F(l,NTR,K,l)* 

VCUSoAM(N)  f XT  R A ( *.  , L ,NTK  , K)  *F  ( 2 , NTR  , n,  l ) *S1  NOAM!  N)  * 

SCOSPhl (N)  ♦ XTkA(3,L,NTk,lO*F(3,NTR,k, L ) *S 1 NGAM ( N ) * 

VSlNPMt(N)  ) 

IF(L  ,NE.  i)  SULN  A ( J , L , N ) = SOLN  A ( J , L ,N ) * S ( NT « , K ) 

GU  TU  710 

3b  SOLN  A ( J , L , N ) = SULN  A(J,L,N) 

S ♦A(M,J,K)*(XTRA(1,L,NTR,K)*F(1,NTk,K,D* 

iCOSGAM(N)  «-XTRAU  ,L  , NTR  , X)  *F  ( 2 , NTR , N , l ) * S 1 NGAM l N)  * 

VCoSPh* (N)  +XTKAI 3 ,L ,NTR,K) *F ( 3, NT  R, K, 1 ) * S 1 NoA.’M  N ) * 

$C  INPH  I < N)  ) 
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$ *Di)FXP(-!N*WAVE'Ic!*S(NTR,K)*XVALlNTF-U  J 
ll(L  .Nfc.  I)  SUL  N -i  ( j iLiM  = 5 OLN  A ( J , L ,N  ) * SI  NTS , k ) 

aO  CuNTlNUE 

r. 

7 ouNTlNUfc 

J.J  SO)  N=1,NRP 
JC  900  L = 1 » 3 

T A = < U. 0.0.0) 

DU  730  J = i, NR  MODE 

1 F ( F.  . f,L  . NTR  ) Go  TO  9 b 

T t>  =CUEXP(-IM*aAVENG*S(M 

in  = ta+SOLN  A! J, L ,Ni *TB*F IL .M  , J , Z i 

iFII  . it.  I ) 7 A = Ta/SIM  ,J) 

UL  TL  73 0 

<0  Ti1  =CDEXP( IN*WAv£No*S (K  , J)*1XVAL(M  ) - KriUjj 
TA  = TA+SCLN  A(J»L»N)*T8*F<L»M  ,J,2) 

I F ( L .NE.  1 ) T A = TA/  S<M  . J ) 

73*  CONTI  NUE 

TA  = TA*CGNST/J:>wRT  IuS1MRHu/lR  4i/l  ) 

TGBL  = TA  vCOEX?  (IM  * *AVE  NO  * KHC  ) 

CALL  MAoANG  (Tool,  TgMAG.  TUANC) 

1 SMAG  = T l)N  AG 
TSANC  - TOanG 

TSDb  = b.bR3890  * OLuG  ( T SMAG  * l.OEo) 

R(1SU8)  = KHO 
UBIL.N.ISUBJ  = TSDB 
ANGU.N.lSUEi  = TSANG 
saved  (M  i = xvaunrslAd  -u 

IF  { M00(ISUo,2)  .Ew.  1)  VllL.NfNn  = U 8 ( L , N , 1 SoB  ) 

1 F I MO  DI I SUC  » 2 ) .Ew.  0)  Y2(L,N,NT)  = DB(L,N,ISU8I 

I M riCOl  ISUB, 2 ) . E W . i ) ANG 11  L » i\ » NT  » = Afw  ( L , N,  I SJ3 ) 

1 F t MOD  ( 1 SUb  , 2 ) . t w . 0)  ANG2(L»N,NT)  - ANG(L,N»ISU8) 

901  CONTINUE 

WHO  = »HD  + DEI  RHG 
iSUb  = ISUb+1 
MP  = M 

IF  (RHu.lE.RHu  MAX)  GO  TO  600 
ISUb  = ISUb- 1 
I F I NT  . NF . NTMAX)  RETURN 
OG  93 0 N= l.NRP 

PGAMMA  = GAMMAIN)/* .7933..9D-2 
PPHl  = PHI IN)/1 .79b3i 90-2 
PRINT  910 

PRINT  927, PuAMMA, PPHl .ZT.ZR 
PRINT  iZU 
PRINT  923 

*27  FORMATl*  GAMM  A ( OEG  ) = • t Ft)  . 1 , • PH  I l OE  G I = • * F o.  i , • <-T  ( KM ) = • , F 10 . 3 , 
%•  Z R ( KM  I = * » F 1 C . 3 ) 

QU  930  JJ=  I iNTMAX 

PRINT  908 » SAVE 0( JJ ) , Y 1( i ,N, JJ)  , ANG1 1 1,N,JJ),YIIJ.N,JJ)  , 
i ANG1( 3 ,N, JJ ). Y l ( 2,N, J J) f ANGi (2,N.JJ),Y2(i,N,JJ). ANGt  l It N, JJ) • 
$ Y2(3.N,JJ).  ANG2 I b . N , J J ) »Y2(2»N,JJ) * AnG £i2,N.JJ) 

90b  FORMAT!  • * ,F7.i»o(IX,P8.3,IX,F7.2)) 
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930  CONTINUE 

If  ( 1 l'PLL'T  ( i » .fcu.u)  RETURN 
CALL  HC PI  TS 
-iETUR  N 

9 i o ConNAf I 1 i 1 i 

i • ELECT:-. if  FIELD  b I (v  FiV  UT  i-t  AS  A FUNCTION  uP  T K A NSH I T T -T  ER  M I N A T jn 

* ul  ST  ANU  ( U i * ,//) 

920  FORMAT  ( * * . 2o  X , • F I t Ll,  AT  RHC.M  I N » , 34  X , * F I E LU  AT  RHuMAX'J 

92  S FORMAT  ( ' ' * 3X  » ' J ' »iX»2(2X|'c2(DB)  • » 2X  , *E2(  ANG  ) ' , 2X , • c Y l Jd ) • .* X , 

* ' FY  ( 4No)  ' ,2X, • tX<  Dd ) • »2X»  ' Ex(  ANU)  * ) ) 

t NO 


M 


1 


SbokubT  IuE  MCFL02 

illicit  heal 

C COMPUTE  t-  I h L S F R OM  n hL  MIN  TO  kHU  MAX 

^ AT  Ufcl  fs  ML  iNTtkVALS. 

c 

CLMVi  M/T^hM/NTtNTk 

o0Mk[\/MLl*\PT/THFTA1^5,pJ*F0Fk(iijf5)txTKAIa»3»23»&l»TUkHlI1iDi, 
i )n/ALI^l  ,FhhQ,KnOMAX,«hJMIN,  JcLkrtJt  JELTAX,cPSLUM2bJ  , 

t.  S I GM  A I 2 3 > ,.MkSLA6tNKMuOfc'»NTMAX 

L UN  ML  * / NC  J 1 Ij  K / A ( Z S » 3 f 3 ) 9 SI  23  9 3)  y c ( 2 3 t 3 ) 9NThS*J(23)  » K VR  A jT  9 l\  V k A T T 9 
i AVRXl  *T  9hV>KTT  iCD^STfUMEGAtWAVti^U 

uuMMLn/KC  K CT/K  {‘♦00)  9 J3 1 j 9 ^ 9 4 J J J , A.MGI  jf  t,  40'3J  , I JPLQT{  l J J , I b JB » NP  P 
tuMMLN/Hl  GN/F  4 J iZ).  t.9^  < 

CCMMON/hb  I NP  T / u AMM  A ( *♦ ) ,PHim,ZT,ZR 
1 9 S INi(«A  l ‘t)  9Co5i>A^l<ti  9 S I i\Prt  I l tf ) 9'-USPHl(  ♦) 
vLAL**i  R.CE.AMo 
v-  OM  PL  tl  X * i o SOLNAI  d,3,  h) 

^UMPlEX*  lo  ThFT  Af  At  S9  C 9 x TR  A 9 Td9TL)dl9  1M/(0.J09  1.D0)/9TA9 

1 FOFR.F ,NTnSC 

■J  t A l <••  b K V r>  A U T 9 K.  V i\  A T T 
J AT  A tr A0/o.37j!33/ 


PRINT  910 
I SUE  = 1 
ubMAX  =■  -1CC0.0 
kHO  - RFC  I'1 1 'J 
,1  = N TR 
X = k HO  - i . 0 
030  CONTINUE 

f Of'  Ik  (FHu.lt.X)  bU  Tu  720 
JO  710  N=i,NRP 
DO  710  L=  1 , 3 
DU  7i0  J = I91NKMODC 
SuLN  A ( j 9 L ,N)  = 10. O9  0.0) 

DO  710  K.  = 1 9 NR.100E 

i F ( to  .i.E.  NTR)  GU  Tu  3 a 

S GL  N A ( o 9 L 9 N)  = S CL  N A(J,L,N) 

i ♦•A(M9J9N.)*(XTkA(  IfLtNTKfKJ’XFlifiNTk  9^91)* 

l CCS  GAM  IN)  +XTF  Alt  ,L  91NTR9  K)  *F  ( 2 , NT*  , K,  i)  *SINGAto(NJ  * 

fcbUSPFl iNJ  ♦XTRA139I  ,NTk,M*F(3»NTk,K,l)'*SINGAM(N)* 
iSTNPHA ( N) ) 

I F ( L *NE.  i)  SULN  h(j  9L9U)  ^ GULN  A(JrLiN)*3iNTR«i\j 
bQ  TL  7 1 u 

3a  SOLM  AIJ,L9N)  = SCL  N A(J,L,N) 

S •*-A(M9J,K)5«ttXTRA(l9L9NTK,K)*Ftl9NTkfK9l)* 

*CUS GAMIN)  i-XTkAU  9L  9 NT R , K ) *F ( 2 9 NT R , K , 1 > *S1  NGA.-U  N)  * 

>CUS PH  i ( N ) ♦XTRA(3,L  ,NTR  ,K)  +F  ( 3 9 NT  k , K,  i)  *SINbAM(N)  * 
iSINPHl (N) ) 

» *CL>tXP(  -lM*rtAVEN(J*S(  NTR  , K ) *X  VAL  ( NTk-  1 ) ) 

IF(L  .N£.  1)  SQLN  A(J,L,N)  = SULN  A ( J , L tN ) * S ( NT  0 , X ) 

I 1 0 u CN T I \UF 
H = fi  — l 
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A = 1 • Jt  O 

II-  (h.GT.O)  X = XVALIM1 
GO  TU  7 Cn 

Uk,  GUN  1 1 NUt 

Utl  N=1,NRP 

l)u  900  L - l » i 

T A = (0.0,0.01 

i)ij  / Ju  J = l,ixKMuOfc 

1 K ( n:-*-  i . ImF . Ni T -* ) GU  TO  46 

tp  -coexpt  -i  m*» a veno*s (mu  , j i*r hui 

1A  = TA+SOLN  A( J , L ,NJ *T b*F ( L ,M+l , J» 2) 

1F(  L .NE.  I l T A ^ r A/  SIMU  , J ) 

GO  TP  730 

9b  Tb  -COEXM  I M*V»AVENU*S  (M*i,  J)  *(  XVAL  l 1)  - PMOl) 

TA  = T A>SGL  N A(  J,L,N)*T8*F(L,M*-i,J,2> 

1KI  .NE.  li  TA  = TA/S(M  + i,J) 

730  CONTINUE 

TA  = TA*CCNST/DSQAT  (uSINIRllG/ErtAO)  i 
TL)bL  - TA  *CDtXr’  ( IM  $ *,7>\/t  Nu  * kHGJ 
GALL  MAGANG  ( TObL  » TDMAG,  TDANGi 
TSMAG  = TOMAG 
TSANC  = Tl)ANG 

TSOo  = d. 605890  * PLUG  ( TSMAG  * l.OEo) 

K ( I su*  ) = P HO 
03 ( L * N * I SOB  ) = TSD3 
A\G ( L ,N,  I SUB  I = T S A No 
90y  G0NT1NUE 

k H(  I = KHU  + UEL  RHO 
I SUb  = I S U b ♦ l 

IF  I F HO  .LE.  BHOMAX*  GU  TO  600 
i SUb  = ISUB-I 
OU  930  K= 1 » NRP 

P&mMMA  = G«MMA(N)/i.7  953290-  2 
,VHI  = PHI  <N)/1.  74:>3^9D-2 
PR  I nT  927,  P GAMMA  , PPtu  ,ZT,ZR 
PP I NT  92  6 
PkINT  929 

927  PUP  MAT ( »C' , LiX , *GAMMa( UEG) =• ,Fd . 1 , ' PM 1 10EG J = » , F u . 1 , • ZT(KM)=*, 

$ P i 0 • 3 » ' ZR(kM)  = ' , FA  0.3) 

928  FURMAT(22X,'Z'  ,37X,'X*  ,37X,' V ) 

929  FORMAT ( CX ,3 (9X , 'RHO(nM) • ,3X, • AMP ( 06) • , 3a, * ANG( Gbo) ' ) ) 

JU  930  L * I SUB 

)>P1  ,\T  9 0tt  ,K  ( J J , OB  ( 1 ,N  , J I , ANG  (L,  N,  J ) , R ( J 1 , 00  l 2 » N , J ) , ANG I 2 , N , J ) * k ( J I 
s , OQ  t 3 , N, J ) , ANG ( 3 * N * J i 

930  CONTINUE 

tF  < IDPLUT  ( 1)  .EQ.O)  RETURN 
CALL  MCPLT2 
kET UP  N 

uOo  rQRMATI  2X,3  l 7X , T 10 .2 , F10  .5  , P 10. '+ > ) 

9 i 0 FORMAT ( • 1 • , 10X, 

i.  'FEgCTkIL  FIFLO  S1KENG1H  AS  A FUNCTION  OF  RHG',//J 

tND 
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Subhi  sj r ine  mCplT  s 

l Mcf'LTb  OrUEx  ATFS  TWO  PLuTS  iFIELO  AMPLITUDE  IN  DtJ  ABOVE  A 

HliruJ  vtLT  Pc*  METER  f ok  i i\W  RADIATED  PURER  VtRSUS  TRANSMITTER- 
TE  »'  i 1 N aT  U \ DISTANCE  FOK  Two  RECEIVER  POSITIONS). 

Cilf-Ht  .*/ T h KM/NT  i NT  F 

wC.MiMl  './SPLOl/S»*VfcUi‘t^)*YL(3,tt,Aj£)  » Y 2 l j » *t  * ‘♦Oa  ) 

o OM  Mi,  iv/  F _ Pl  OT/R  IAOO)  , OBI  jfHt  •♦JO)  , AN  G I 3*  •* , 40  T ) > I DPLdT  ( 1 J ) * ISJd*NKP 
CDMwr  N/X°LUT/XMiN,  X INC  , YMlN,  V INC  > S I ZEX,  SI  ZEY 
bUi  ML  N/  mG  1 NP  T / ijAMM  A ( A ) ,PHI(4),ZT  , Z R 
l , SING  AM(4)  , CCS GAM  I A)  , S INPhI L ) , LOSPH I ( 4) 

COMML  N/MClNPT/ThETA(c5»:>),PUPk(23»:>)»XTRA(J,3»23,5)»TJPrtT(<::>)» 

1,  XVAL12  5)  , FKEw»M)UMAX,RhGMINf  DEl*HU , OELTA X, EPSl  ON  1 25)  , 

i,  S I uM  A l 2 3 ) , NR  SC  A 6, Nk  mOD  E * NT  MAX 

C CM  PL  t X*  l o THt  T A , F OFK  t X Tk  A 

K L A L < b XVAI,FRPG,RHLMAa*RHGMIN,OELRHG»DElTAa»EPSlUN,SIuMA,TGPHT 
*<  r A t X=  j S i NG AM , C 0 SG  A M»  SINPH1  ,COSPHI 
RLAl*6  GAMMA, PHI  ,zt  ,ar 

REAL  XLURVF (21/0. , l./ , YCURVi(2)/2*0./,YCJKV2(2) /c*.2/  , 
i YL0RV3i2)/2*.4/, YC0RV4L2 )/2*.o/ 
j I MENS! UN  Y(40U) 
l>  ! MFr.S  I ON  GAMMAIX  A ) , Phil)  (A) 

CoF'PLtXf  i«  LUMP  (3) /•*.  COMPONENT  • , *X  CUMPuNENT 

<.  ' Y COMPONENT  V 

00  900  1-1,2 
Ou  900  IL'EGIN  =1,3 
v-ALL  uGNPL  t IBFGIN  ) 

GUI  Y AXmNG  ( ^ . ) 

CALL  I NT  AX  j 

CALL  PAGE  (SI2FX  + 3.,SIZEY  + 3.) 

LAL  L PHY. 'OF  ( 1 . 3 i L . 2 ) 

CALL  TITLH'  i,  'TRANSMITTER-TERMINATOR  L)ISTAnCE(NM)  ' , J3  * 

» ' UL  ABOVE  I Uv/M  FUR  l KW  * ,24, S12EX,SI2EY ) 

C*»L  L GRAPH!  XMiN,  X INC,  YM  IN,  Y INC) 

NO  cOJ  J=1,NRP 

oAMMA  t)(  J)  = GAMMA(J)/1. 7453290-2 


>’HJ  0(  J)  = 

PHI L J ) / 1 

. 7‘*33«_9D-2 

IF  ( J 

• L*  W • 

1 ) 

LALL 

RESET l • DASH* ) 

if  ( J 

• by. 

2 ) 

CALL 

DOT 

I F ( j 

• by* 

o ) 

CALL 

CHNUOT 

1 F ( J 

• b*Q. 

4) 

CALL 

DASH 

DO  5C0  R= 1 , NT  MAX 

L F ( I .EG.  I)  Y(K)=Yi(  IBEGIN. J,n) 

I F ( 1 .EG.  2)  Y( K)=Yc( IBEGIN, J,K) 

3 00  CON  T 1 NOE 

call  curve ( savfq.y , nt max,o ) 
on  Continue 

CALL  MESSAGL IUPLUT,4C,1.  ,9.3) 

CALL  MESSAGLCONPt  IbEGIN),l6, L . ,9 . I) 

CALL  MtSSAGCFFEw  = * ,7.»1.,<J.9) 

CALL  kEAlNOL FREG, 3, 1.7, fa. 9) 

CALL  MES5AGC2T  = ZR  = • , 17 , l . , to  .7 ) 

CALL  RFAI NU(ZT,2,i. 7b,8. 7) 

LALl.  w t ALNC  ( ZF  ,2,3.40,6.7) 
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1 


91'' 


CALL  ME  SS  AC ( • k£CE  iVfcR  DISTANCE  = *,20, 1.0,0.31 
iH  I .tu.  i 1 CALL  KfcAL  vOUHUMINtl  (3.3,8.31 
I L ( l « Fl.  . *_  1 CAL  L k F AL NG(  RHUMAX  , 1,3.5,0.51 
CALL  «tES^  Ao  ( • C AMMA  = PHI  = • . 20  . t . 3 . 9 . L 1 

CALL  * E AL  M,{  GAMM  AO  l 1 ) » i , 7 .0  » 9. 1 1 
CALL  fv  L AL  Nu  ( PH  1 ')  ( i I f l , o . 3,  * • X 1 

o«L  L ME  SS  AC  (•  GAMMA  — PHl  - 1 m.0  i o.D  i o .9) 

CALL  K E.AL  NO(  GAMMAD1  £ 1 fit  7.0,  0.91 
C ALL  kEAL  M<IPHiD(2  1 ,1 ,0.3,  o.91 


LALL  ML.  SSACt  • u AMMA  = Pul  = • t c 0 t o . u t o . 7 1 

CALL  l’. F.AL Mj  ( GAM MAO ( 31 , 1 , 7.0,  o.  /l 
CALI  REAL  N(.|PhlL)(3l  , i . o.3,b.7) 

lALL  ME  S S Au  ( * G AMM  A = PHI  - 'iiUiL.0iO.il 

o AL  L kE  AL  NC  i GAM  MAC  ( A 1 » L » 7 . 3 * d . 1)  I 
L.ALL  KtALNO(PHlUU)  tl  ,8.3,8.51 
CALL  ENDGM1) 

CALL  UkEL ( 9. 7 t o . bo  1 

CALL  TiTLEC  'fit'  ' » U,  ' '(Jfi.fi.) 

CALL  Ck A P H ( 0 • ii.iO.ti  .1 
CA.L  L HE  Sr  T I 'DASH*  ) 

CALL  CukVE  ( XCU»  VE  , YCoRV<t,2  ,01 
CALL  DOT 

CALL  CUKVEIXCURVE ,YCURV3,2,C1 
CALL  CHNDOT 

CALl  C UPVFIXC UR VE,YCOkV2»2, 01 
lalL  JASti 

CAL  L CURVE  UCURVE,YC0RVLf2t91 

CALL  LiVCPL  i leEuIN  I 

CUNT  1 NJE 

RETURN 

tND 


I 
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iUbk  uT IKt  MCPLT2 

L *CPcT2  GENERATES  ONE  PLoT  (HELD  AMPLITjOC  IiN  UB  ABOVE  A MICRO  VOLT 
C (>F  ' 'IE  T t R FOR  1 RAoluTtU  FOi.cn  WbkSoS  UISTaICE  FROM  TkANbMI  TTcK)  • 
v.  H Ml  :*  / MC  PL  CT  / R ItOO)  , i)0<  3 ,t.40J)  , A N G ( 3,4,40  0)  , I DPL01  i 131 , 1SUC,NRP 
■ : Ml  './  XPLl  T/XM  In  , X In  C f YM1  \ , y I nC  , S I ZE  X * Si  Ze  Y 

C U M ML  >\! / H G I N P T / G AM M A ( *♦  ) * P n I ( *t  I • 4 T , Z R 
t , S 1 KG  AM  I 4 i » CUbGAM ( 4)  , S INPHl  <41  , ..USPHI  ( 4i 
C Gi'IMl.'J/  MC  INPT/THETa(4.3»3)  ,EUFR(23,o),XTRA(;»,3,23,3)  , TUPHT12  3)  , 
j.  XV  AL  ( 2 3)  , FRL  J ,RHUMAA  , KriUM  IN,£JELRhC)*UELTAX,EPSlON(23)  , 

F S 1 GM  A ( 20 1 *NkSLAB»  NKrtUOE  * NT  MAX 

COMPLEX *10  THtT  A,  FOLK  ,XTKA 

4 c A L * d AVAL  » FREw,  KHOHAX  , RriiJK  IN,  OElRHG  , DECT  AX,  £P  St  IN,  S I oMA  , TUPhT 
^FAL*b  oAMMA,  Ph  l , ZT  ,ZR 
* t A L*  6 S1NGAM.CUSGAM,  SlNPHUCOSPhl 

REAL  XcURVf  (21/  3.  , 1 . / . YC JRV1( 21/2*0. /,YoJRV2 <21  /2*.2/ , 
i YCUKV312) /2*.4/,YCU RV4(2)/2*.o/ 

DIMENSION  Y ( 400 1 
DIMENSION  GAMMAH4)  *PHI  D(  41 

cOMPLEX*U  COMPlil/'t  COMPONENT  •<  COMPONENT 

i *3  COMPONENT  '/ 

GO  9 JO  lHEolN  =1*3 
-ALL  boNPL ( IocGIN  > 
v-ALL  YAXANCO.l 
CALL  INTAXS 
CALL  PAot (ll.,xl.) 

CALL  T IT  L L ( ' • , -l,  'KriO(KM)  •*  7,  *0b  AfeOVE  t UV/M  FuR  i k.*',<l4, 

* SIZLX.S1ZFY) 

CALL  uPAPh(  AM  In,  XINL,  YMIi\,Y  INC) 


5uO 

oO  J 


GO  oO 0 J= l , NR  P 

v,AMmaC(J)  = GAMMA ( J ) / 1.74a  32  9 0-2 
PHluMJl  = FHI  ( J 1 / 1 . 7*,  332  90-2 
1F(J  .EU.  11  CALL  RE  S ET  ( ' DASH* 1 
I F ( J .Eg.  2)  cALL  DOT 
I F { J .Eg.  31  CALL  CHNDOT 
1 F ( j .FQ.  4)  CALL  13  Ab  H 
JO  pCC  K= 1 , I SUo 
Y ( K J = DB<  IBEGIN, J,X) 

CONTINUE 

CALL  CURVE(R*Y* 1 SUB , 0 ) 

CONTINUE 

CALL  MESSAG( iCPLJT ,40 ,4 . ,9.3 1 

CALL  MESSAGI CUMP(  IB EG  INI, lo, 4. ,9. 31 

CALL  MESS AG( • FREQ  = ',7, 4. ,9. 11 

lALL  k E AL  NCi(  F RE  0 » 3 » 4 . 7 * 9. 1 1 

CALL  NtSSAGCZT  = ZR  = • , I /,  4.  , 8.  91 

CALI  REALNtHZT  ,2,4.84,8.9) 

CALL  KEALNO(Zk,t, 0.34.  ,8.9) 

c AL  L MESSAG( • GAMMA  = PHI  = *,20, 4. ,0.7) 

CALL  REAL NO(oAMMAO( l 1 , 1 , b .0, 8 . 7 1 
CALL  REALNU(PHID( 1 ) ,1 ,6.3,8.71 

CALL  -IESSAGI  • GAMMA  = PHI  = ',20, 4. ,8. 3) 

c A 1 L kcMLNUI GAMMA0(2j  ,1,3.0,0.31 
CALL  * E AL MU (Ph 10(2) ,1 ,0.3,3.31 
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CAL  L 

MtSSAGL  • GAMHA  = 

PHI  = 

4 v 2 J f H • t t>  • 3 ) 

v,  AL  l 

REAL  NG(  CAIMADL  JJ 

v i i 5*  Ot  8«i ) 

CALL 

,<  E AL  NO  ( PH  I Q(  3)  li 

9 0 • f dm  JJ 

CAL  L 

'IcSSAGL  ' GAMMA  = 

PHi  = 

v-  A L L 

R E «l  Nu(  u AM ‘4  AD  ( i 

9 i f !>iOf  & • i ) 

^ L L utALNutPHlL/i*,}  | 1 f t)  «bf  d • Li 

CALL  ENO&L'.ill 

CAL  L PtiYSOR  (d. ,9. ) 

CALL  T I T L t 4 • ',1.'  SO,'  * . D , L . , l . i 
>, AL  L orvAPnlO*  * 1 • i 0 « » 1 • ) 

CALL  kbSFT ( ' 9 A S H ' ) 

v-ALL  CuRVEi aCUR VE,  YCoRV<»,2»01 

v.ALL  JOT 

CALL  CU*VL(XCURV£,YCURV3,2,0) 

CALL  CHNUOT 

CAL  L CURVEL XCURVE,YCURV2.2.0) 
o al  l dash 

CALL  CJRVE(XCURVE,YCuRVi,2»0) 

CAL  L tNl'PL  ( IdLoIN  ) 
r i nu  f 
return 

END 


t 
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iUbPLOT  INF!  MCSTEPCM) 

C CJM'*UT  r jtUti>  aL  I it  C ML'Uc  CONVERSION  CO  EF  F LC  I EN  T S . 

IMPLICIT  REAL  *rf(  A -n,  U-Z  ) 

CC-Pl  N/TEXM/NT  * NT  K 

CiH^i  \|/CAP/CAPI  (Xb,b,  ■>)  . TM  JXM  ( X P « P * 5>J 

LijMM  ./mC  1MPT/THET  A ( p , p ) , FORK  t 2 P * p ) , XTkAlj*i*2:J*  p I * T DPrlT  l <.  0 ) # 

» XVAL(XbI*Fk£U  *k  hJi-IAa.  , F hL)MlN*JELkHU*OELTAx*r.  P S L U 1 M le-'Ji  * 

i S I GM A ( 2p ) ,NRSLAb,NRMODE,NTMAX 

Cl#. 'A  Ml  ./  ACS  TOR  / A ( 2 p 1 p,  p j • if  2p  * o J iC  < 2 t»  * b ) * N TH  Su( 2 P ) *X VK  A.)l  * KV  AATT  > 

* AVKK,UT»».V0-  isTT»CO'JSr*iJi'/ILOA»wAVCi'#iJ 
C J'i  fV  *_■  X*  1 1»  C D t X P 

C Ur!  r’L  E X*  l o nTIiSw 

COM  PL  EX*  it#  THFTA  , FiJFk  , A , S,C  , TNUaM.CAP  I , 

* l,M/(  O.OO*  L .00)/  ,8f  b)  ,ANSl  p)  , TS I 5 , b ) , XTRA 
k t:  A l * ■ # K V f-  A 0 T ,nVRATT 

i\  E A L * N F R K 


-IP  = M+  l 

Ic^  . lC.  NTR)  RETD  Pi* 

00  1 7 N = 1 * NKMUOL 
17  bl M ) = ( 0 . * 0 • i 

1HMF  .hC.  NTh)  GO  Tu  Xi 
jl: ; «.y  k = 1 , nf  mode 

Oil  ii  L •-  1 * Nh’^lJUt 

00  ->0  J - 1*Nki1UDE 
TSfL.Jj  = TN'}PM(M,  <_♦.*> 

Oj  p(L)  = bill  *A(  MP  * J *n  ) *COtXP  I -1  •'1*hA  VE  N0*$  I MP  * J I MXVALIM)  ~ 
L XVAL ( M P ) ) ) *C  AP  l ( M,  L*  J ) 

CALL  CLINFu(TS»p*  AN  S * NR  MuO  t *p*o*E-<K  ) 

00  c.7  f = i.NPMUOE 

27  A(M(i*K>  = ANSI  )/S(NTK,n) 

Du  it  N=1,NKMUI)E 
lb  rtlN ) = (0  • • 0 • ) 

2 9 COMF1  N'JE 
GO  TP  X-* 

XI  DO  Xj  K = 1 , NX.-lGDt 
00  ».b  L L=  l * NRMc  Ub 
DO  Zb  L = l.NRMJOt 
T 5 ( L L * L I = T NORM! M * L L * L ) 
cp  o(L ) = CAP  1 ( M * L * K ) 

CAlc  Cl  I N EDITS* b*ANS»NRMODE*5*0*ERR) 

DU  33  J = 1 * NR MODE 

ip  A ( 0 * J , X I = ANS ( J ) * S ( M * J ) / S 1 N T R * K ) 

2 J CONTINUE 
2*»  CONTINUE 
RETURN 
END 
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SUBKoUTlNE  MOHNKL  ( Z , H 1 , HZ , H 1 PRME , H2PRME I 
COMPUTE  MUU1F1EJ  rlAUKEL  FUNCTiUNS  UP  URDEH  ONE  THIRD 
IMPLICIT  RIAL  *8  ( A-h  ,0-Z  ) 

COMPLEX*!!}  COS  JRT  » CDt  XP 
R t AL  * j COAoS 

COM  PL  EX* l o Z,  I , HI , HZ  , HI PRME , nZPArte , ZPUw ER , T ERM 1 , T ERMZ  » 

% T=KM3,ZTLkM,  TERM,  SUM  1 , SUMZ  , SUM3  » SU.M*+  »SGKTZr>* 

f EXP  a ,l:XP4,  Ex  P3,  EX  P4.EXP5  ,GM2F,GPMFP,MP0wCR,  iE  T A,  R TZ i 

$ C0NST1,CC.\ST2  ,COJ  ST3 , C0NST4 

01MENSIUN  AU3li  ri(2j)t  C(23l,  0(43),  CAP(i*t) 


U 1 ME  NS  I UN  A (43)  , 8(23), 
OaT  A A/ 

l.  9. 304867io9300C0D-0l 

i>  3. 7434jo3242 34300  02 

k 5.4l6b5437404340D  02 

> 2.oo2o33l8707‘tOOD  01 

) l. 840127394410000-01 

!■  2. 913341420000000-04 

i 1.  41557C0CC0J000D-C7 

!>  2.O0000000000000D-11 

OAT  A 6/ 

» 6.  782987251400000-01 

l i.  19629404787J500  02 

» 7. 474221821372000  01 

; 2. d5325737403000D  00 

> 1. 0422939453 000 00-02 

. 2 . 25262E90CO0  3000-03 

I.  9.  7290C0000  00  5000- 09 

> 4 . JOOOOOOOOJOJOOD- 12 

jat a c/ 

> 4. o54l83  58‘,o0000D-0  1 

. 5. 2213059311 40 OOD  01 

. 2.  706427187022000  OA 

> 9.  ibl 3000*5109000-01 

> 4. 84244l03b0000C0-03 

> o. 19659900000000D-0O 

i 2. p^bOOUOOOOOOOOD- 09 

i 0. COOL OC09GG090CD  CO 

OAT  A D/ 

> fa. 7d<.98725140000D-01 

i 1 . I9o244047b73300  03 

t 1 . 420 1 02 1 4609d6 50  03 

» 7. 989l20fa47290000  01 

» fa. 0 76487783230000-01 

i 1.036301278000000-03 

> 3.33C74COOCOO0OOO-O7 

. 1. 02 CuO 00 00000 OOD- 10 

DATA  CAP/ 

> 1 . 04 1 666666600 6 7D- 0 1 

» 2.91 b 49 U 264641 400“ 01 

> 1. 49937fa29bb862far  01 

i 3. 20  74  9CC910000CO  03 

I.  1.  79190200000000D  Ob 


3.101‘*5:> 723097000  01, 
b.70217o3pl900800  02, 
2. 379454463 oiJ200  02, 

0.  1210,'0439C3oC00  00, 
2. *6 3303 096400000-02, 
2.562749500000000-03, 
8.870000000000000-09, 
1.0000 00000000000-12/ 

1 . 1 3 0-+ 9 7 o 73*40000  Ol, 
1.533710317786500  02, 
3.235393862 152000  01, 
o* 136 03 73635100 CD-Oi , 
2. 105503122000000-03 , 
l.9l3o7lOCOCOOOQO-Co, 
3. 890000000000000-10, 
0.000000000000000  00/ 

6.20291 1446190000  00, 
6.215840394215000  01, 
1. 121301940796000  01, 
1. 9 i 28 A 2 63 43 900 0D-01, 
o.O5oo3o82OCO3COQ-04, 
3. 1oj300u00000000-07, 
l . 300000000000000-10 , 
O.OOOOOOOOOCOOC 00  00/ 

4. 521941500 96 20 uO  01, 
1.993823413144500  03, 
7.118306496735100  02, 
1. 902171582688000  Ol, 
8. 42402  0489600000-02 , 
9. 3 do 78o 900000000-05, 
3.4l 3500000000000-08, 
5. 00 CO 000000000 00 -12/ 

8.3530347 22222220-02, 
8.8162  726  7443  7580-01  , 
7.892301301 13870D  01, 
2 . 40 6fa3 960000 000  04, 
1.746437700000000  07/ 


2.067o37  A-»d7Jl6C0 
6.287787l442bo400 
9. 345b495Joo3i 000 
a.  l392803o4*t3G00D 
2. 8842  Ob  01 0000000- 
2 • 02  3o8o OUOOJC 000- 
3. OIOOOJOOCUOOOOO- 


5.3833234A34ja000 

1.27809l93l486d00 

1.078531287384000 

1.0937o7b00980000- 

2.331677830000000- 

1.444700000000000- 

3. 200000  0 JO 0 0000 o- 


4. 58434ot359l5000 
4. 8 7 3 1 o 8 9 300  39000 
3 . 394357  3 j2  330000 
3. 3l2229oo9400  000- 
6.355018200  00000  0- 
3. 322000000000000- 
8. 000000000000000* 


3.766326230801500 
2. 044947090382 060 
2.69o32b2iS4u0i0D 
3. 71 8 3 10 523339 OCO 
1.002o2l4S690u000- 
7. 512435000000000- 
1.96200000000GOCD- 


1.282203  7435o3270- 
3. 321408281862770 
4. 74451 33od680000 
1 • 989231 2wO00C OOD 
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DATA  l/('.DCtl.D0i/ 

0 >li  K 00  f 3/  i • 7 3 A OAU  Ok-  1'jut  ttbD  JO/ 

D*.T  i.  ALf’HA/ti.i>3oO/Aloo3H'S:>l0-Ul/ 

_»  AT  ^ 'w  • ’is  S T 1 / t i?.jiic;l'»0<tPl0<:>22iJ-01  »— S.o3,<43i2oAJ'»Jk>7u-Jii/ 
jaI  , CiNS1«./(  t » •* . c.  j ii.  j 3Ao2  3*7  Jo  70“  J l ) / 

0 A T A CONS  Ti/(-'#.jbv2jo2o2J^Oo7c(-Jlt  A.3odlvJ-*:>i,JA'jAAG-.»iJ/ 
.> AT  lUNSTh/  < — j'itjCiutd'jyo/J-uii'i.aoJH  Jt3l  J2aJ20-^i.  )/ 


ZPGwf  K-l  .0 
bL«Ma=  o.c 
SuM4  = J . 0 
ZMAG=CQAFSm 

1F(ZMAG  .GT.  4.2J  GG  TO  70 
I F { AN  AO  . C»E . J.A)  GO  TO  10 
M=U 

ou  Tl  3C 

1F(  AM  AG  . Gt  . *♦ . I ) GO  TO  20 
N=  1 3 

GO  T L J G 

\ = 2a 
SuMl=0. 

3UM2-0. 

tTEKM=-A**3/200.J 

Ob  IQ  M=  1 , \ 

aUM  i=  Sui».  1+A(M)*4PGaF* 

SOM  A=  $ u M2  ♦ B ( M 1 * Z P Ok.  E F 
a uM  j = SO M 3 ♦ C ( M ) * L P U r>  C * 

S UM 4-  SUK4«-D(M)*  aPCaE* 

4Pr>FA  = ZFL'FER*ZT£RK 
1 F ( C L AES ( Z POWER ) .Lfc.  1.00-301  Gl 
GuN  1 1 NOE 

G M2 F = 1* ( Z* SUM A- A. *SUMl ) /ROOT  3 
UP,VFF  =I*(SUM4*2.*Z*Z*SUM3)  /ROOT  3 
rll=Z*SUM2+GM2F 
ri  c — h 1 — <•  • 0 ^ G 1*1 2 F 
HiPf  M£=SUM4+GPMFP 
H?P«<ME=F;IPRmF-2  .O'FGFmF P 

r.F  T Ur< 


TO  oO 


SOM  1=1.0 

sum  4= 1 . c 

l<TZ=CDSWKT  l Z) 

S0PTZd=RTZ*Z 

ZTFPM=I/S0RTZB 

mpcufk=i.o 

T E R M = — 1 . 3 / Z 
DU  B(J  M-=l»l4 
AP0aER=ZP0WE3*ZTPRM 
MPL  aE K-MPGVsER* (-ZT EKM ) 
f\  RFi  1-CAP  ( M ) * Z P 0 W E R 
T Fp Mi  =C  AP 1 M ) EMPOWER 
S Uvi  * - S U M l T £ R M 1 
SUK^=  SuMcMERMA 
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i.U<4'  = SUMj  + K*TFPMl 
i U '4  *«  = S U M 4 ♦ M * r E K 4 2 
60  CONTINUE 

aUtoi  = iUN3*TC*f-! 

i uM  <*  - b 'J  M *+  * T fc  K h 

EaP1  = cDEaHZ.*I  *SQk  J <.0/3.1 

cXP,.=  EXPl*C  ON ST  1 

EXP^CuNiTc/tXPl 

EXP4=CJN3T3*EXP  1 

bXPii=GON3T4/EXPl 

l> ETA=  ALPHA/ CDSURT  (RTZ  ) 

<.RE AL=Z 
Z I MAO=- 1 *Z 

If  IZREAL.Gfc  .O.I).Ok.<.IMAG.GE.O.O)GU  TO  90 
rtl=BETA*( EXP2*SUM2*FaP5*SUM1  * 

HiPhME=tiETA*(EXP2*(SuM2*(-0.2:>/Z'*-I*KTZ)'*-SUM4>+EXP!>*(SU/41*(-0.<::>/Z 
$ - I *R  TZ ) + SUM3  ) i 

Gli  TG  11C 

90  Hl=iitTA*EXP2*SUM2 

HlPkME=6tTA*EAP2*lSUM2*(-0.2f>/Z*-I*KlZ  I ♦SUM  4) 

110  IH  (Zr<EAL.GE.G.O.OR.ZlMAG.4.T.O.O)GO  TO  120 
H2=  BE T A* ( EXP3*SUM1+EXP4*SUM2) 

H2P  *<ME=  BET  A*  ( EXP3*<  3uMl*(  -0 . 2 !>/ Z - l *RT  Z ) +SUM  3 ) +E  XP  **  ( SUM2* l -0. 23/Z 
1 + i *RTZ ) ♦ S UM4 I) 

aE  T UP  N 

UO  U2»Lf  TA*fcXPi*Svl-4l 

U2EKMt=BcTA*EXP3*lSJ(4  1*1-0. 2a/Z-I*RTZ  ) ♦SUMj  i 

RETURN 

END 
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SUdhCUT  IN  E MAGANG(  At* G, MAG,  ANGLE  I 
I ►.P  LICIT  tv  E AL  *8(A-h,0-Z) 
kIAl^cj  USUI'.T,GAklL'S 
o b M RL  E X * i t>  ARG,  IM/(  0.  DO, 1 .00  )/ 

K c A L ^ a MA  a 

JATa  ROTl>tG/57.Z9:>77S  Si/ 

£\T  hV  f tGNGL  F lARG,MAG»  AwCl  E ) 

AKGRAL  = ARG 
At-GMAb  = 1M*Ak0 

.AAG  - L- Sw f*T  ( Ak GK A L A AR  GR  A L * AkGMAG* AKGM AG ) 
IFlK.Ab  .EL.  3.0)  Gb  1U  10 
t,bSL  — AkGRAL/MA*, 

IE  (LUSG  .lT.  -1.0. AnD.  COSO  .GT  . —1.01  ) CUSG 
IFtCbSU  • uT  . 1.0. AND.  COSO  .LT.  1.011  CUSu  - 1 
L>  ANGLE  = DAR  COS  ( CPS*  ) ♦■ROT  jtu 

iRtAkbMAG  .LT.  0.0)  hNGLE  = 3o0.0  - ANGLE 
RETURN 

10  COSO  - 3.0 
>0  TC  5 
uND 
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dc'RALuriivfc  nrcAlMZ) 

U C GMPUT  fc  fci  » = a » r V Hfc  1 uHT  GAINS  FOR  TRAN  SHI  T TtK  ANO  RECEIVE. 

IMPLIlIT  COMPI  EX*  lo  (A-H.O-A  I 
C GHMLiV  rloT  l.  MP/F  F l ( c 'j  » o ) * FFi  ( 23  » j > 
lUMMCN/hTGN/F i ^ ,2  St  t>,  2 ) 

tOHMCJ.N/MC  INPT/THETAC2  j,  b j , FCFk  l 2 , 3 J , XTPA  l 3 , 3 , 2 5 , 3)  » T LiPH  T 1 2 5)  . 

1 X V AL  ( 2 3J  » F REG  »K  HGMAX  * R hOM  I N , jE  l RnU  fOfcLTAX*tPSLliN(23)  . 

i SI GMA ( 2> ) t NKSL  A 3,  iSIKMGDt  , NTHAX 

CuHMl  N/MCSTUR/  U23»5,3J  , S 1 25 , 5 i , C ( 25,  3J  ,NTHSG ( 2 5 1 ,K tfR  AuT,  KVKATT  , 
ft  AVKROT  »AVRKTT  tCONST  » UM  fcG  A , wAVfc  .kQ 

CUMPlfcX*Lo  GDSwKT 
R fc  AL  * 3 DfcXP 

COMPLEX*  l <j  NGSw  . I M / (O.UO»  i .0  0)  / ,NTHSw 

is  E A l * d XVAl.FREu.KHOMAX.RHLjMI.M,  OELR HU » OfcLT AX . fcP SLON , S I GMa 
R E A L*  8 KVP AOT,KVP ATT  » AVRkJT , AVRKT T , CONS T , GMfc ^A , W AVENU 
P E AL*  8 L ( 2 ) , EPSON  0/6.  B3A3  AO- 12/ , ALP  HA/ 3 . I aO -*t/ , F AC  l 
KpAL*d  KSuP 

REALMS  TOPhT 
DO  10  D M~  1 1 NR  SL  A3 

NGSw  = IfcPSLOMM)-lM*SlGHA(M)/QMEuA)/EPSLNO 

DO  100  K=  1 » NRMUOfc 

SSO  = S(M,M**2 

SQROOT  = CDSOKT(NGSw-SSO) 

cSU  = ClM»N)**2 

RSwR  = S.ROGT 

IF(RSs(K  .LT.  c.l  S0RgU7  = -SQKUUT 
00  100  I Zs  1 1 2 

>j  = kvkatt*  icsy-*-ALPHA*z  ( iz) ) 

= KVRAT7*CSJ 

CALL  MDHNKL ( 00 , Hi  0 , H*. 0, H 1 PRMO , H2PR M 0 ) 

CALL  MDHNKLi Q,Hl,H2fHlPRM,H2PRM) 

k,AP  HI  0 = HI  PRMO  + AVRKT  T*H10 

CAP  H^O  = HZPrMO+A VRKT  T *H20 

FAC  2 = I M*K VRAQT*SuRuOT 

PAC3  = PAC2/NGSw 

FI  = -tCAPH23-FAC3*H^01 

F2  = CAPHIO-F AC3*H1 0 

F 3 = - ( Hii  PF  MO-F  AC2*H«_  0) 

FA  = HIPKMJ-F AC2*H10 
F AC  1 = DEXP( ALPHA/2. *Z( III J 
FlhMfKdi)  = FAC1*(F  l*HUF2*H2) 

F ( 2 »M  , K * I Z ) = ALPHA/I  I M*2. ** AV ENU ) * F ( l,M,K,  IZ  )♦  l . / l M*  A\/RKUT  *F  AC  l * 

$ (F l*HiPRM+F2*H2PRM) 

F(3»M»K»Ia)  = F3*Hl*FA*H2 
FFI ( M » K ) = F1*H10+F2*H2C 
F F 3 ( M » K ) = F 3*H 1 Q+F  A*H2  0 
Fll,M,K,IZl  = F(1,M,K, IZ)/FFl(M,K) 

F I 2 » M f K , I Z I = F{2,M,n, IZ»/FF1(M,k) 

P ( 3 » M » K » I Z ) = F13 ,M,n, IZ)*FOFR(M,K» /FFj(MtM 
100  CONTINUE 
RET UF  N 
6N0 


riE  •** 
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SUbf-CuTlNE  CLIN  EU  (A,  6,  X,  N, 

S N DIM,  IFLAG,  EFF) 

C 

L Lt  IN  USES  L-U  DECCtor CS I T I ON  TC 
„ I 1 ,U  T nr  T ^ ] a NuU  L Ak  M AT  kICES  L*  U 

C.  Such  T 1 1*.  T L * U - A.  C AND  L‘  ARE 

C bljkFu  1 N A.  ThLS  FUkM  lb  USED  wITh 
C oALK-SU^ST  IT  uT  l ON  TU  FIND  THE  SOLN 
L X OF  A * X = L * U * X = 6. 

C f,  IS  Ti-t  NUMBER  OF  EOLATIONS  AND 

C Li  1 M lb  THt  DIMENSION  OF  AlL  ARRAYS 

C IN  THE  P ARAMt TE K LIST. 

C 

C IF  IFLAG  = 0,  L,  U,  ANU  X ARE 
C COMPUTED. 

C IF  IFlAG  IS  NUN-ZERO,  IT  IS  ASSUMED 

c that  l and  u f*ave  been  cgmputed  in 

C A PREVIOUS  CALL  AND  ARt  STILL  STUKED 
C IN  A.  TnUS  ONLY  X IS  lUMPU TED. 

C E kk  lb  THE  ESTIMATED  RELATIVE 
C ERrsUk  l.f  THE  SOLUTION  vECTJK. 

c 

CUM PL  EX* lo  A,  8,  X,  T 
! NT  EGEK*2  I ROW 
JIMENSICN  AIN  DIM,  N DIM), 
i BIN  U 1 h J , XIN  DIM) 

DIME! SIGN  IRUW(BC),  JlBC) 

DA  i A cPS  /i.Oc-l'u/ 

0 

c 

IF  (N.GT.oO)  GU  TO  vDO 


IF  1 1 FL  AG .NE 

.0) 

GU  TU  600 

DO  OB''  I = 

0(1)  = 0.0 

l , N 

DO  040  J = 

l , N 

GG  = oUAdS  (All, 

J)  ) 

J4Q 

IF  U(  I)  .1 T. 

vvl 

Gil)  = GU 

IF  It  I I ) .EG. 

J.JJ 

ol  TO  GO  I 

DBO 

oUNTJ  \Ut 
iPk  = cPS 
PPIV  = 0.0 
DO  iOO  I - 

1 ,N 

100 

1 R G W ( I ) = I 

OH  BOO  L = 

1 » N 

PIVOT  = 0.0 
x = L - 1 

DO  240  I = 

L , N 

IF  IK.LT.l) 

GU 

TU  2 j 0 

DO  22  C J = 

1 .K 

220 

All  ,L)  = All 

,L) 

- A ( j , L i * A(  1, J) 

z JO  *-  = C DAB  S ( A ( I , L ) ) / wll  ) 
IF  IPIVUT.GT.F)  GU  TU  2*0 
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f 


P l V (j  T = E 
NPIVUT  = I 
2 4^  CONTINJF 

IP  (PlvGT.tc.0,0)  ( »u  TG  901 
It-  (PPI  V.LE.P  1VUT  ) *0  TO  260 

EkR  = ERR  * PFlV  / PIVOT 
ir  (ERR.GE.l.O)  GO  TO  901 

2 60  Pf  I V = PIVOT 

IP  (NPI VGT.EJ.L I GG  TO  2dO 
UtNPlVOT)  = OIL! 

J = I ROWll ) 

IkOkiL)  = I RUh ( NP I VUT ) 

IRUwl  ivPIVQT  ) = J 
GC  2tO  I = l , N 
T = A( L,  I J 

A<L  , I I = AINPIVGT, I ) 

A J NP I VOT » 1 1 = T 
2t>0  CONTINUE 

260  IF  (L  .E0.NJ  Go  TO  6J0 

T = ( l.GDC,G.GDo>  / h (L  (L) 

K = L ♦ 1 
M = L - 1 
DO  430  I = K,N 
IF  (M.LT.i)  GO  TO  400 
00  iuU  J = l.M 

3 60  AIL, I)  = AIL, I)  - AIl,JJ  * AIJ  ,ii 
400  AIL  ,1  I = T * AIL, I ) 

4 60  CONTINUE 
6 00  CGNI I NUE 

IF  (ERK.GT.l.OE-6)  PRINT  998,  EPR 
C 
C 

600  00  o20  I = 2 , N 
o20  Xll)  = (0. (.00, 0.000  | 

J = 1 KOk  I 1 I 

XII)  = 81 J)  / All,  11 
DC  7C0  I - 2 , N 
J = I ROW  I I ) 

K = I - 1 
DO  t>60  L = 1,K 

660  XlII  = XI  I J + AII,U  * X I L ) 

XII  I = (BIJJ  - XII))  / All, II 
700  CONTINUE 
k = N - l 
DO  61  0 I = 1 *K, 

J = N - I 
M = J ♦ l 
(JO  bCO  L = M , N 
XlJl  = XIJ1  - X ( L I * AIJ, LI 
300  CONTINUE 
kETUFN 
C 

900  PRINT  999 
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err  = 1.0 

x E T Ut-  N 

901  HkIM  997 

e r-  r = i . r 

K f T U r \ 

9 9/  r'JkMM  {'lERKUR  IN  Cl  IN  Ew>  MaTkIX  IS  SINGULAR  J 
rQf  Mil  ( 1 CAJT  lu  \i-'  i 

'>  ' lI  IN  EU  HAS  JECuMRUSEO  AN  I a -CONO I T I JNE  D M4Ti\IX.  »/» 
t,  ' 1-taulTS  "ILL  -iAvt  RELATIVE  ErRO« 
rukMT  CilKWOR  IN  LlIN  EU,  MATRIX  SIZE  uh  EATER  ThAN  jj*  i 

END 
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Mission  Research  Corporation 
Attn:  Dr.  C.  Longinire 
One  Presidio  Avenue 
Santa  Barbara,  C'A  93101 


Pacific-Sierra  Research  Corp 
Attn:  Dr.  E.  Field 
1456  Cloverfield  Blvd. 

Santa  Monica,  CA  90404 


Mitre  Corporation 
Attn:  Mr.  G Harding 

Dr.  John  Morganstern 
Ms.  Susan  Morin 
P.O.Box  208 
Bedford,  MA  01731 

Lockheed  Research  Laboratory 
Attn:  Dr.  J.  B.  Reagan 
Technical  Library 
3251  Hanover  Street 
Palo  Alto,  CA  94304 


(3)  University  of  Denver 

Denver  Research  Institute 
Attn:  Mr.  Herbert  Reno 
Mr.  Donald  Dubbert 
Colorado  Seminary 
University  Park 
Denver,  CO  80210 

(2) 

University  of  California 
Attn:  Dr.  Henry  G.  Booker 
P.O.  Box  109 
La  Jolla,  CA  92038 


Los  Alamos  Scientific  Laboratory 
Attn:  Dr.  J.  H.  Coon 

Dr.  R.  F.  Tascheck 
Dr.  D.  R.  Westervelt 
P.O.  Box  1663 
Los  Alamos,  NM  87544 

Boeing  Co. 

Attn:  Mr.  J.  F.  Kenney 
Mr.  H.  R.  Willard 
P.O.  Box  3707 
Seattle,  WA  98108 

EG&G,  Co. 

P.O.  Box  384 
Kenmore  Station 
Boston,  MA  02215 

President,  Sandia  Laboratories 
Attn:  Mr.  W.  C.  Myra 
Albuquerque,  NM  87115 

President,  Sandia  Laboratories 
Attn.  Mr.  T.  B.  Cook 

Dr.  Byron  F.  Murphy 
Livermore,  CA  94550 

Lawrence  Radiation  Laboratory 
Attn:  Dr.  Frederick  D.  Seward 
Dr.  Glenn  C.  Werth 
P.O.  Box  808 
Livermore,  CA  94550 


(3)  Pennsylvania  State  University 

Attn:  Ionospheric  Research  Lab. 
University  Park,  PA  16802 

University  of  Illinois 
Attn:  Aeronomy  Laboratory 
Dept,  of  Electronics  Engineering 
(2)  Urbana  Campus 

Urbana,  IL  61801 

Geophysical  Institute  of  the 
University  of  Alaska 
Attn:  Technical  Library 
College,  AK  99735 

Stanford  Research  Institute 
Attn:  D.  Nielson 

Dr.  E.  T.  Pierce 
Technical  Library 
Menlo  Park,  CA  94025 

Radio  Science  Laboratory 
Stanford  University 
Attn:  Prof.  R.  A.  Helliwell 
Stanford,  CA  94305 

Rand  Corporation 
,,,  Attn.  Dr.  C.  M.  Crain 
1700  Main  Street 
Santa  Monica,  CA  90406 
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R&D  Associates 

Attn:  Dr.  Forest  Gilmore 

P.O.  Box  9695 

Marina  Del  Rey,  CA  90291 

Science  Applications,  Inc. 
Attn:  Technical  Library 
1 200  Prospect  St. 

La  Jolla,  CA  92037 


